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On the basis of an analysis of the numerical results corresponding to the half-filled one-dimensional 
t-t'-V model on some finite lattices, Kokalj and Prelovsek (KP) have in a recent paper [Phys. Rev. 
B 78, 153103 (2008)] concluded that the Luttinger theorem does not apply for the Mott-Hubbard 
(MH) insulating phase of this model (corresponding to the region V ^ t) in the thermodynamic 
limit; KP even suggested, incorrectly, that failure of the Luttinger theorem were apparent for a 
half-filled finite system consisting of A'' = 26 lattice sites. By employing a simple model for the 
self-energy E oi a. MH state, we show that the finite-size-scaling approach of the type utilised 
by KP is not reliable for the system sizes considered by KP: although the Luttinger theorem is 
exactly satisfied for our model, by employing the latter finite-size-scaling approach and the same 
lattice sizes as employed by KP, we obtain comparable quantitative amounts of "violation" of the 
Luttinger theorem as reported by KP for various values of V/t in the range [3.6, 10]. On the basis 
of the equivalence of the model under consideration (at half-filling and for small values of t' /t — 
ideally, for t' — 0) and the XXZ spin-chain Hamiltonian for SU(2) spins, we further show that 
for V greater than a critical value Vc{t,t'), where Vc{t,0) — 2t, the system under consideration 
has a charge-density-wave (CDW) ground state m the thermodynamic limit, corresponding to a 
doubling of the unit cell in comparison with that specific to the underlying lattice. Although 
this ground state is also insulating, its spectral gap is due to the broken translational symmetry; 
it is not a correlation-induced MH gap. The Luttinger theorem (as generalised by Luttinger for 
inhomogeneous ground states) is therefore a priori valid for this broken-symmetry ground state. 
This fact rigorously establishes that the conclusion by KP is indeed erroneous. Finally, we present 
a heuristic argument due to Volovik that sheds light on the mechanism underlying the robustness 
of the Luttinger theorem. In an extensive appendix, we present the details of the calculation of the 
single-particle Green function of the broken-symmetry ground state of the model under consideration 
(for i Vc) by means of bosonization and in terms of the form factors of a class of soliton-generating 
non-local fields pertaining to the quantum sine-Gordon Hamiltonian, as determined by Lukyanov 
and Zamolodchikov (LZ); the latter fields can be shown to represent the slowly- varying chiral fields 
of our bosonized problem. 

PACS numbers: 71.10.-w, 71.10.Pm, 71.27.4-a 
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I. INTRODUCTION 

Considering for simplicity a system of spin-less elec- 
trons, for the uniform iVc-particle ground state (GS) of 
this system one defines the Luttinger number 

iV, = ^e(G(fc;/.)), (1) 

k 

where 0{x) is the unit-step function, G{k;e) the zero- 
temperature limit of the thermal single-particle Green 
function in the grand-canonical ensemble in which the 
mean value of the number of particles is equal to N^, 
and /i the zero-temperature limit of the associated chem- 
ical potential (see later). Thus G(fc; e) corresponds to the 
A^o-particle GS of the Hamiltonian H of the system. The 
summation with respect to fe in Eq. (1) is over the entire 
available wave- vector space; for systems defined on a Bra- 
vais lattice, the latter region consists of the corresponding 
first Brillouin zone (IBZ). The Luttinger theorem under 
consideration states that [1, 2, 3, 4, 5, 6, 7, 8]: 

N^^N,. (2) 

There are some aspects to be taken into account (in par- 
ticular in dealing with finite systems) concerning the def- 
inition of 0{x), specifically at x = 0, which have been 
elaborated on in Ref. [7] and which we shall not discuss 
here. 

Inspecting the proof of the Luttinger theorem [7] , one 
concludes that the validity of this theorem is dependent 
on the equality of the value of the chemical potential 
fi on the right-hand side (RHS) Eq. (1) with that of the 
thermodynamic variable p, specific to the grand-canonical 
ensemble to which G{k;e) corresponds. With Em-o de- 
noting the energy of the M-particle GS of H, one readily 
verifies that the implicit dependence of G{k; e) on the 
thermodynamic variable fl is solely through the require- 
ments 

En^-o — En^-i-o < p, < En^+i-o — En^-o. (3) 

For metallic A'o-particlc GSs, where is necessarily 
macroscopically large, the deviation of E^^^i-q — E^^-q 
from Epf^-o — Epf^-i-o is microscopically small, of the or- 
der of (see Appendix B in Ref. [7]). For these states, 
the value of p, is therefore up to an infinitesimal correc- 
tion uniquely determined and is equal to Ej^^^i-q — Ej^^.q. 
The latter quantity coincides, up to an infinitesimal cor- 
rection, with the zero-temperature limit of the chemical 



potential fi/j, where /3 is the inverse temperature, satisfy- 
ing the equation of state corresponding to as the mean 
value of the number of particles in the grand-canonical 
ensemble under consideration. Consequently, neglecting 
infinitesimal corrections, for metallic GSs the value of 
the fj, on the RHS of Eq. (1) is bound to be equal to the 
zero-temperature limit of iip. 

For insulating A'c-particle GSs, En^j^-i-o — En^.q devi- 
ates from Ejsf^-Q — E^^^i-o by a finite amount, equal to 
the magnitude of the fundamental spectral gap. Since 
for these states, and for p satisfying the inequalities in 
Eq. (3), the value of G{k; e), Vfc, e, is independent of the 
precise value of p, one naturally concludes that the Lut- 
tinger theorem in Eq. (2) should apply for all values of 
fi satisfying the same inequalities as p. That this is not 
necessarily the case, was first discovered by Rosch [9] in 
his study of a simple model of a MH insulator. 

In Ref. [7] (see also Ref. [8]) it is unequivocally demon- 
strated that the possibility of the breakdown of the Lut- 
tinger theorem for MH insulators has its root in the pos- 
sibility of the existence of a false [10] zero-temperature 
limit of a sum over Matsubara frequencies, and that this 
false limit is eliminated through identifying the fi on the 
RHS of Eq. (1) with the zero-temperature limit of /i/j. 
Thus, insofar as the /x on the RHS of Eq. (1) is concerned, 
insulating A^o-particle GSs are to be dealt with in exactly 
the same way as their metallic counterparts (see the pre- 
vious paragraph). In this connection, it is important to 
realise that from the standpoint of the Lehmann repre- 
sentation of the thermal single-particle Green function 
^# (of which G is the zero-temperature limit), metallic 
and insulating systems are fundamentally similar at non- 
vanishing temperatures [7]. 

We should emphasise that were it not for the possibil- 
ity of the existence of the above-mentioned false limit, 
the Luttinger theorem would apply for all values of fi 
satisfying the same inequalities as p. Neither Stanescu, 
Phillips and Choy [11] nor KP [12] appear to have appre- 
ciated this fundamental aspect. 

II. THE t-t'-V MODEL IN ONE SPACE 
DIMENSION 

The system considered by KP [12] is defined on a regu- 
lar one-dimensional lattice and is described by the Hamil- 
tonian 

N 

H = Y,[-itc]cj+i+t' c]cj+2) + ll.c. + Vnjh,+i], (4) 
where 

Uj = cjcj (5) 

is the occupation- number operator of site j. In Eq. (4), 
the site with index j ~ N + j' is identified with that with 
index j = j', where j' = 1,2. Unless we indicate other- 
wise, in what follows we shall use the lattice constant a 
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as the unit of length so that the set of k (hereafter k) 
points relevant to the sum in Eq. (1) is a subset of the 
interval (— tt, tt] (the points — tt and tt are identified along 
the k axis). Using the periodic boundary condition, for 
the lattice under consideration, and with N even, one 
has: 



IBZ = 




For the half-filled case, considered by KP [12], 

No = In. (7) 

A. Discussion of the numerical results by KP 
concerning finite systems 

The finite systems for which KP [12] have carried out 
explicit numerical calculations, correspond to 

iV= 14, 18,22,26,30, (8) 

accommodating respectively Nc = 7, 9, 11, 13, 15 spin-less 
electrons. For these systems, KP [12] employed exact di- 
agonalization of H, making use of the Lanczos algorithm. 
We shall later show that in the thermodynamic limit, and 
for V > Vc{t,t'), the half-filled GS of the H in Eq. (4) 
is not uniform, but a non-uniform CDW state [13] (see 
Sec. Ill herein). This broken-symmetry GS cannot be 
reached by extrapolating the finite-size numerical results, 
which by necessity all correspond to uniform GSs. 

In Fig. 2 of Rcf. [12], t{l+V/t) Re[G(fc; ^i)] is presented 
as function of the discrete values of fc, namely k = 2ttI/N, 
I = 0, 1, ... , N/2 (in Ref. [12], fi is chosen as the origin 
of the energy axis). We note that G(fc; /z) G R, Vfc [7], so 
that the explicit use of Re[G(fc; /x)] may imply that in the 
calculations reported in Ref. [12], Im[G(fc; /i)] ^ 0. This 
is relevant, since the numerical artifact Im[G(fc; /i)] ^ 

can in principle give rise to a false violation of the 
Luttinger theorem [7] (see in particular Sec. 6.2 herein). 

From the data in Fig. 2 of Ref. [12] one observes 
that for iV = 26 the Luttinger theorem, Eq. (2), in- 
deed applies for all V/t e [0,8]: at fc = 2Trl/N, for 

1 = —6, —5, . . . , 0, 1, . . . , 6 (that is, for 13 out of 26 dis- 
tinct k points of which the underlying IBZ consists), one 
has G{k;fi) > (cf. Eq. (1)). Notwithstanding this fun- 
damental fact, KP remarked that [12]: 

"we note that /cl is indeed near 7r/2, however, 
even without finite-size scaling a small devi- 
ation /cl 7^ tt/2 may be observed for V > At J' 

Here is the Luttinger wave- vector, defined by KP [12], 
and others, as the solution of 

G{k^;^i)=0. (9) 

In Ref. [7] (see in particular Sec. 2.4 herein) it has been 
emphasised that this and similar equations for higher 



space dimensions (defining the so-called 'Luttinger sur- 
face') must be used with caution, for these equations 
may not have solutions. In the particular case at hand, 
where the underlying IBZ consists of a finite number of 
points, and where specifically tt/2 ^ IBZ for iV = 26 (cf. 
Eq. (6)), it should be abundantly clear that "deviation 
/cl ^ 7r/2" is devoid of any meaning, whatever. It seems, 
rather surprisingly, that the criterion adopted by KP [12] 
for the validity of the Luttinger theorem amounts to the 
requirement of Eq. (9) having the exact solution k^^ = it /2 
(and, by symmetry, = — 7r/2), even in the cases where 
7r/2 does not belong to the underlying IBZ! We repeat 
that, for the reasons spelled out in Ref. [7], even in the 
limit N ^ oo, for which the IBZ consists of the contin- 
uum set (— 7r,7r], absence of a solution to Eq. (9) cannot 
be held as signifying breakdown of the Luttinger theo- 
rem. 

One general technical remark is in order. In Rcf. [12] 
KP indicated that 

"it has been realised that the LSR [Luttinger 
sum rule] should as well apply to finite sys- 
tems [10,11]", 

where "[10]" refers to Ref. [14] and "[11]" to Ref. [15]. In 
this connection, we emphasise that this result is already 
implicit in the original paper by Luttinger and Ward [1]; 
nothing in the proof by Luttinger and Ward is dependent 
on whether the sums with respect to fc in the underlying 
expressions run over a finite or a macroscopically-large 
number of k points, to be approximated by a continuum 
set; the only significant aspect associated with these sums 
is that they include the complete set of the k points that 
are relevant to the uniform GS under consideration — 
otherwise, by using an incomplete set of k points, the 
invariance property of the trace operator under the per- 
mutation of the product of the operators on which it 
operates, which is vital to the proof of the Luttinger- 
Ward identity, is violated [7]. What in our opinion has 
contributed to the obscuring of this very basic fact with 
regard to the applicability of the Luttinger theorem to 
the uniform GSs of finite systems, is the explicit reliance 
by Luttinger and Ward [1] on a continuum model; in 
contrast to lattice models (more explicitly, those with- 
out physical boundaries), the uniformity of the GS of the 
continuum model adopted by Luttinger and Ward is inti- 
mately connected with this model being infinitely large. 

Summarising, on the basis of the data as presented in 
Fig. 2 of Ref. [12] we conclude that: for at least N = 26 
and V/t G [0,8], the explicit calculations performed by 
KP [12] undisputedly confirm the validity of the Luttinger 
theorem, Eq. (2), specifically in the MH insulating phase. 
If this observation is applicable to all values of N pre- 
sented in Eq. (8) above, it follows that the main conclu- 
sion arrived at by KP in Ref. [12], that the Luttinger the- 
orem broke down in the MH insulating phase, must en- 
tirely rest on the outcome of the finite-sizc-scaling scheme 
adopted by the authors. Following this observation, be- 
low we critically examine the last-mentioned scheme. 



B. An analysis of the finite-size-scaling scheme 
employed by KP 

In order to examine the reliability of the finite-size- 
scaling method adopted by KP [12], we employ the fol- 
lowing model for the self-energy S(k; e) at e = fi: 

S{k;ii)=A+-^^^-^, ke{-TT,Tr], (10) 

where A and B arc constants to be determined. Evi- 
dently, by definition A and B arc functions of /i, however 
for our following considerations knowledge of the func- 
tional forms of these parameters will not be necessary. 
We emphasise that the expression in Eq. (10) is hypo- 
thetical and is not based on a microscopic calculation. 

With Ek denoting the non-interacting single-particle 
energy dispersion, for which we choose the tight-binding 
expression 

Efc = -2<cos(fc), (11) 
for G{k; fi) one has (in the units where h ^ 1): 

G{k-ii)^—^, (12) 

where 

Ek=ek + S[k■^J). (13) 

We remark that use of a more general expression for the 
energy dispersion e^, one involving t' in addition to t, 
would only increase the number of parameters in our con- 
siderations without affecting the main conclusion of this 
section. 

By employing a graphical presentation of , one can 
determine a continuum of values for A, B, t and fi for 
which the GS corresponding to the Green function in 
Eq. (12) is a half-filled MH insulating GS. In Fig. 1 we 
present the Ek and the corresponding G(fc; ^) for the 
specific values A = 1, B = 1.33, t = 1 and /i = 0. 
Note that our specific choice /i = is consistent with 
the explicit particlc-holc asymmetry in the calculations 
by KP [12], where t' ^ (sec the diagram in the upper 
panel of Fig. 1). 

Before proceeding with our analysis, two remarks are 
in order. Firstly, for the reasons indicated in Ref. [12], 
for the system considered by KP, G{k; /i) must be small, 
Vfc e IBZ, for large values ofV/t, scaling like t/V^. This 
is clearly not the case for the function G(fc; /i) displayed 
in Fig. 1. One can directly convince oneself that this as- 
pect is of no consequence to our conclusions to be arrived 
at in this section. To this end, one should consider the 
fact that k = ±7r/2 remain the exact zeros of G{k; ji) 
on replacing E^ and /Lt by aEk and a//, where a is a 
non- vanishing but otherwise arbitrary positive constant. 
This implies that for in particular ^ = 0, through a sim- 
ple re-scaling of the parameters A^ B and t, one will be 
able to render G{k\ ji) as small as desired for all k £ IBZ. 
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FIG. 1: (Colour online) The correlated single-particle energy 
dispersion E^, Eq. (13), and the associated G{k;^), Eq. (12), 
both in the units where a = 1 and h = 1. The parameters 
used are, in the latter units, t — 1, A = 1.6, B — 1.33 and 
fj, = 0. From the behaviour of Ek one observes that for a finite 
range of values of jj, indeed the underlying GS is insulating; 
this is the case for /i = (broken horizontal line). In the lower 
panel, the small circles in the vicinity of fc/vr = 1/2 mark the 
values of G{k; /i) at fc = 7r/2 ± n/N for the values of A'^ listed 
in Eq. (8). Although for the model considered here, /cl/tt is 
exactly equal to 1/2, the 'finite-size scaling' scheme based on 
the last-mentioned values of A'^ yields /cl/tt = 0.500742. See 
Table I. Note that indeed G{k; fj.) > for fc £ (-7r/2,7r/2) 
(cf. Eq. (1)). 

One can readily verify that such uniform scaling of the 
above-mentioned parameters cannot lead to any change 
in the conclusions of our analysis to be presented be- 
low (with the aid of the details that we present in ap- 
pendix A, the reader may verify this statement through 
explicit calculations). One should in addition note that 
since the smallest N that we consider in our following 
investigations is equal to 14, Eq. (8), we never probe 
G(fc; /i) farther from k = 7r/2 than 7r/14 « 0.22 on ei- 
ther side of this point. Consequently, for our specific 
considerations in this paper it is entirely irrelevant what 
non- vanishing magnitude G(fc; /x) has outside the interval 
[7r/2 — 7r/14, 7r/2 -I- 7r/14] (and by symmetry, also outside 
the interval [-7r/2 - 7r/14, -7r/2 + 7r/14]). 

Secondly, concentrating on the half-filled GSs (see 
Eq. (7)), we assume that the expression for S{k]^) in 
Eq. (10) applies for all N, so that by increasing the value 
of N one merely samples the same function at a more 
dense set of k points of which the IBZ in Eq. (6) con- 
sists. Denoting the IBZ corresponding to a system con- 
sisting of M lattice sites by IBZ(M), with k e IBZ(M) 
and k e IBZ(M'), M ^ M', in reality the value of 
S{k]ii) (and thus of G{k;fj,)) pertaining to the M-site 
system is different from that pertaining to the M'-site 
system. Since in the t-t'-V model direct hoppings take 
place between nearest and next-nearest neighbours, in 
the absence of long-range order one may assume that 
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for M,M' > Mc, where Mc is a finite integer, the last- 
mentioned difference is insignificantly small. The quan- 
titative similarity between the results that we deduce be- 
low and those reported by KP in Ref. [12], leads us to 
believe that the values of N in Eq. (8) satisfy the inequal- 
ity N > Mc. 

Following KP [12], we define 

GiN) = I (G(| + ^; ^) + G(| - ^; /i)) , (14) 



and 



AG{N) - G(| + ^; ^) - G(| - ^; m). (15) 



For the smallest and the largest values of N given in 
Eq. (8), one has respectively n/N sa 0.224 and n/N « 
0.105 (sec the small circles in the lower panel of Fig. 1). 
Using the scaling relationships [12] 



ZiG(iV) = ^ + 4, 



(16) 



(17) 



KP [12] obtained: 
dG{k-n) 



dk 



k=ir/2 



lim 



AG{N) 
2ti/N 



By employing the explicitly-calculated values of G{N) 
and AG{N), corresponding to the values of N in Eq. (8), 
KP [12] deduced the values of the parameters {ai, bi, ci} 
and {62, C2} with the aid of a least-squares fitting method. 

Taylor expanding G(ki,; fj) (for N = 00) to linear order 
at fc = 7r/2, imposing the condition in Eq. (9), KP [12] 
obtained the approximate expression: 



G{n|2■^i) 



k =-- 

' 2 dG{k;fi)/dk\k=^;2 



Clearly, the accuracy of this expression diminishes for 
increasing values of [/cl — 7r/2|. On the basis of the obser- 
vation that ai/62 ^ for V/t e [3.6, 10], KP [12] arrived 
at the conclusion that for the latter range of values of 
V/t, the Luttinger theorem broke down. See Fig, 5 in 
Ref. [12]. 

By employing the expressions in Eqs. (16) and (17) and 
the Green function in Eq. (12) as determined in terms of 
the model sclf-encrgy in Eq. (10) and the Sk in Eq. (11), 
on the basis of the same approach as KP [12] (specifically, 
by making use of the values of TV given in Eq. (8)), we 
have calculated the parameters {ai,6i,ci} and {&2,C2} 
(see appendix A for some relevant technical details). We 
present these parameters corresponding to ^ = 0, < = 1, 
B = 1.33 and some values of A in Table I. In this ta- 
ble we also present the values for ki^/n as calculated on 
the basis of the approximate expression in Eq. (19). One 



TABLE I: The parameters {ai,a2,as} and {02,62} (in the 
units where a = 1 and h = 1), Eqs. (16) and (17), for various 
values of the parameter A, Eq. (10), obtained by means of a 
least-squares fitting approach described in appendix A. In all 
cases, jj, = 0, t — 1 and B — 1.33, in the above-mentioned 
units. For fitting we have used exact data corresponding to 
the values of A*' listed in Eq. (8). In the last column we present 
ki^/n, with as calculated according to the last expression in 
Eq. (19). The values for /cl/tt should be compared with the 
data for the same quantity displayed in Fig. 5 of Ref. [12]. 



A 


ffli 


bi 


Cl 


62 


C2 




1.6 


0.0046 


-0.4255 


91.080 


-12.509 


-99.290 


0.5007 


1.7 


0.0455 


-2.6752 


129.57 


-11.087 


-145.89 


0.5082 


1.8 


0.1145 


-6.3715 


187.00 


-9.0193 


-209.58 


0.5254 


1.87 


0.1900 


-10.340 


244.79 


-6.9890 


-269.15 


0.5544 


1.9 


0.2325 


-12.544 


275.86 


-5.9145 


-299.84 


0.5786 



observes that for increasing values of A, the deviation 
of ki^/ir from the expected value 1/2 increases. It is seen 
that in our model calculations, the increase in the value of 
A brings about a similar change in /cl as does the increase 
in the value of V/t in the calculations by KP [12]. Since 
for the model under consideration the zero of G(fc; fi) at 
k = tt/2 does not depend on the value of A, it follows 
that the deviation of the calculated values for k^/ir from 
1/2 is entirely bogus and a direct consequence of not em- 
ploying sufficiently large values of N. This undisputedly 
demonstrates that the conclusion arrived at by KP [12] is 
devoid of any significance as regards the validity or oth- 
erwise of the Luttinger theorem. In contrast, if KP [12] 
had not disputed the validity of the Luttinger theorem 
and imposed the condition G(7r/2; /i) = as the result to 
be expected for N 00, they would have gained some 
valuable information regarding the behaviours of G(fc; fi) 
and S{k; fx) for k in the close vicinity of fc = tt/2. For the 
reasons presented in Sec. 2.4 of Ref. [7], we should em- 
phasise that satisfaction of the condition G{tt/2;^) — 
is not a prerequisite for the validity of the Luttinger the- 
orem at half-filling. Nonetheless, as we shall see later 
(see, e.g., Eq. (C242)), G(7r/2;/x) = is exactly satis- 
fied at half-filling and under the conditions for which the 
bosonization scheme is exact for the model under consid- 
eration. 



III. A HEURISTIC ARGUMENT IN SUPPORT 
OF THE LUTTINGER THEOREM 

The reasoning that we present in this section is due to 
Volovik [16] and is related here with his kind permission. 

According to Volovik [16], a metal and an insulator, as 
characterised by respectively a pole and a zero in their 
corresponding single-particle Green functions, belong to 
the same topological universality class, or have the same 
topological invariant in the momentum space (Chap. 2 
in ReL [17]). That is to say. the transformation of the 
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last-mentioned pole into the last-mentioned zero at the 
metal-insulator transition does not change the topology 
of the Green function [17]. This implies that the adi- 
abatic transition between the two states, as well as the 
adiabatic transformation between non-interacting and in- 
teracting fcrmions (Sec. 8.1.6 in Rcf. [18]), cannot give 
rise to 'spectral flow': the state without excitation, the 
vacuum state, transforms into a state which is equally 
free of excitations. Consequently, as in the case of non- 
interacting fermions, for interacting fermions the number 
of particles in the GS remains to be determined by the 
'volume' of the Fermi sea in the case of metallic states 
and by that of the Luttinger sea [7] in the case of insu- 
lating states. 

For completeness, topological stability is defined 
(Sec. 8.1.6 in Rcf. [18]) as signifying the property by 
which any continuous change in a system leaves the sys- 
tem in the same topological universality class (see the 
previous paragraph). A continuous change in the sys- 
tem may be brought about by an adiabatic change in the 
coupling constant of the interaction potential. 



IV. THE CDW GROUND STATE OF THE t-t'-V 
MODEL IN ONE SPACE DIMENSION, AND THE 
LUTTINGER THEOREM 

In the thermodynamic limit, the half- filled GS of the 
one-dimensional t-t'-V model is a broken-symmetry state 
for V > Vc{t,t') (to be precise, there are two such GSs, 
differing only by a translation over an inter-site distance). 
FoTVc{t, t') one specifically has T4(t, 0) = 2t (appendix C, 
in particular Eq. (C98)) [13]. With (O) denoting the GS 
expectation value of O per site, for the half-filled broken- 
symmetry GS one has (see Eq. (C63) and the subsequent 
remarks) 



{h,) = - + (-1)-' 5n, 



(20) 



where 5n is the CDW order parameter. In the uniform 
GS, 5n = 0. Thus, the primitive unit cell of the CDW 
GS consists of two sites, leading to the IBZ corresponding 
to this state, which we denote by IBZ', Eq. (B7), hav- 
ing half the size of that corresponding to the half-filled 
uniform GS, Eq. (Bll). The CDW GS is an insulator 
whose gaps in the single-particle excitation spectrum is 
located at the boundaries of the IBZ', exactly as in the 
case of ordinary semiconductors. In contrast to the lat- 
ter systems, whose discrete translation symmetry is due 
to a periodic ionic potential (an external potential), the 
discrete translation symmetry of the CDW state under 
consideration (a Z2 symmetry) is due to a spontaneous 
symmetry breaking (see the second paragraph of Sec. C 4 
on page 22). This breaking of the translation symmetry 
of the underlying Hamiltonian by the GS being possible 
only in the thermodynamic limit, for any finite realisation 
of the one-dimensional t-t'-V model the corresponding 
GS is a uniform one. The half-filled GSs that KP [12] 



have determined numerically are therefore all uniform, 
characterised by 5n = 0. 

The Luttinger theorem specific to inhomogcneous GSs 
is explicitly discussed by Luttinger in Sec. V of Ref. [2], 
under the heading "The Band Case" . The role played by 
G{k\ fi) in the case of uniform GSs (see Eq. (1)), is in "the 
band case" played by the eigenvalues of the Hermitian 
single-particle Green matrix G(fc;/^), which is the exact 
representation of the single-particle Green operator G{e) 
at e = /i in the momentum space. For the CDW state 
under discussion. (G(fc; /^) is a 2 x 2 matrix (appendix B). 
For this state, the statement of the Luttinger theorem 
is that in Eq. (2), in which the Luttinger number A^l is 
defined according to (cf. Eq. (1) above; see Eq. (95) in 

Ref. [2] and note that y^{k] /i) = (/i — Qk^ifJ.)) ) 



(21) 



where {y^{k;^)} arc the (real) eigenvalues of the Green 
matrix G(fc;/i). Owing to the insulating nature of the 
CDW state, in the light of Eq. (7) one immediately ob- 
serves that at half-filling the Luttinger theorem is triv- 
ially satisfied for this state. To appreciate this fact, one 
should realise that for all k £ IBZ', y±(fc;/i) ^ and 
that, since each unit cell for the CDW state under con- 
sideration consists of two lattice sites, for this state the 
number of k points of which the IBZ', Eq. (B7), consists 
is half that of the points of which the IBZ of the uniform 
GS consists. 



V. THE INHOMOGENEOUS GROUND STATES 
OF THE t-t'-V MODEL IN ONE SPACE 
DIMENSION 

In this section we explicitly consider the broken- 
symmetry GS of the t-t'-V model in two limits, one cor- 
responding to F I Vc{t,t') (for t' sufficiently small in 
comparison with t), and the other to V ^ t t'. In 
appendix C we present the details underlying our calcu- 
lations corresponding to the former limit. 



A. The case o{ V > K(t, t') 

For the spectral gap £q corresponding to V larger than 
and sufficiently close to Vc{t,t'), one has (Eqs. (C190), 
(C104), (C226), (C212) and (C214)): 



£c^2M ^ 4V2t cxp 



W2{V/V - 1) 



(22) 



where Ai is the mass gap, Eq. (35). Owing to this ex- 
ponentially diminishing £c, as | V^, for sufficiently 
close to Vc one can employ a continuum description for 
the calculation of the correlation functions of the sys- 
tem under considerations at distances that are large in 
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comparison with the lattice constant a of the underlying 
lattice. With x = ja, in the framework of this descrip- 
tion one expresses the lattice annihilation operator Cj as 
(Eq. (C5)): 

~ ^^ { e+''^'^(^°) R{x) + e'^'^^^J") L{x)], (23) 

where the slowly-varying chiral field operators R{x) and 
L{x) can be represented in terms of the bosonic field op- 



where U± are Klein factors [19] (see Eq. (3.14) herein) 
which ensure that Cj, as expressed according to Eq. (23), 
appropriately removes one particle from any eigenstate 
of the number operator, in the relevant Fock space, on 
which it operates. As we shall indicate in appendix C, 
the signs '±' with which the RHS of Eq. (25) is pre- 
multiplicd (which equally could have been chosen to be 
'=F') is of significance for the preservation of causality in 
calculating the single-particle Green function. At half- 
filling, where fcp = 7r/(2a), Eq. (C3), the phase factors 
g±fc,Oa) ^23) reduce to {±iy . We point out that 

the symbol '~' in Eq. (25) signifies two facts: first, that 
the expression in Eq. (25) applies for a 0, or \x/a\ ^ 1, 
and; second, that it is a leading-order term of an infinite 
Fourier-type series the formal form of which is presented 
in Eq. (14.30) of Ref. [20]. Further, the in Eq. (23) 
becomes a strict equality provided that the Fourier spec- 
trum of Cj, Vj, Eq. (C4), is limited to two intervals of 
width 2yl centred at ztfep, where A is equal to a fraction 
of Ti/a. 

The dynamics of the slowly-varying chiral field oper- 
ators R{x) and L{x) is governed by the quantum sine- 
Gordon Hamiltonian [21] in Eq. (C76) below, in which 
the parameters g, u and K are fully determined by t, 
t' and V (Eqs. (C77), (G123) and (C116)); provided 
that t' is sufficiently small, at half-filling the low-lying 
excitations are affected only perturbatively by t' (see 
Eqs. (C35), (C41) and the remarks following these equa- 
tions, in particular those in the paragraph following the 
latter equation, on page 16). 

Let K'g and K'^,, Eq. (B6), denote two reciprocal-lattice 
vectors corresponding to the Bravais lattice of which the 
underlying unit cells consist of two lattice points, and 
k e IBZ', Eq. (B7). For the {£,£') element of G(fc;z), 
where £,£' G {0,1}, k e IBZ' and z e C, one has the 
expression in Eq. (B9), in which 

Gjj'(^)=.lim Gjj'(w„). (26) 



erators n{x) and ${x), satisfying (in the thermodynamic 
limit — see Eqs. (C18) and (C19) — and in the units 
where h = 1): 

[S{x),n{x')]_^i5{x-x'), (24) 
as follows (Eqs. (C6) and (C27)): 



(25) 

I 

Here (Eq. (25.14) in Ref. [22]) 

G,- {tu;„) = hm r dr e^-^ (r, 0), (27) 

in which (in the units where h ~ 1) 

(2n+ l)7r 
w„ = ^ , n = 0,±l,..., 

P 

is a fermion Matsubara frequency, with P the inverse tem- 
perature (in the units where fee = 1), and (Eqs. (23.3) 
and (23.6) in Ref. [22]) 

^j.j'ir.r') - -Tr{poT,[c,(-^r)4(-^T')]}, (28) 

where is the grand-canonical statistical operator, 
the r-ordering operator, and 

T = it (29) 

the imaginary time (similarly for r'); Cj{—iT) is the 
imaginary-time counterpart of Cj in the Heisenberg pic- 
ture (similarly for c],(-ir)) (Eq. (24.1) in Ref. [22]). 
In the zero-temperature limit, the trace on the RHS of 
Eq. (28) reduces to the GS expectation value of the r- 
ordered product of Cj(— ir) with c^, (— ir'). We note that 
in Eq. (27) we have identified r' with zero on account 
of the fact that in equilibrium ji (r, r') is a function of 
t-t'. 

In the CDW state, Gj,ji{iujn) is a non-trivial function 
of j and j' , in contrast to the uniform state where it is a 
function oi j — j' . With reference to Eq. (22), for ], 14 
the deviation of Gjj'{iuJn) from a function of j — j' be- 
comes exponentially small so that, for V sufficiently close 
to Vc, to an exponentially small error one can approxi- 
mate Gjj'^iujn) by a function of j — j'. On doing so, and 
by renumbering the site indices according to the scheme 
in Eq. (Bl), in what follows we shall restrict our con- 
siderations to the single-particle Green function G(fc; z). 




exp [- iy^ / dx' n{x') ± iy/^^{x)\ , 
J — oo 



8 



defined according to (appendix B) 
G{k;z)= e-*'=(J") 0^,0(2), where fc e IBZ. 

] = -N/2+l 

(30) 

Following the considerations in appendix B, the function 
G(fc; z) as defined here can be viewed as consisting of the 
union of Gofi{k; z) and G'i^i(fc; z), in which k e IBZ . 

Let A{x, —it) and B{x, —it) be the imaginary-time 
Heisenberg-picture counterparts of the Schrodingier- 
picture operators A{x) and B{x), where A{x) and B{x) 
may be one of L{x) and R{x). Defining, in analogy with 



Eq. (27) (in view of the Gjfi{z) in Eq. (30), below we 
suppress x' which would be the analogue of /), 

Gi,B(a:,ic^„)= lini / dr e'^"^ ^^^^{x,t), (31) 

where (cf. Eq. (28)) 

^A,B(a^,T) = -TT{p,,Tr[A{x,-tT)BHO,0)]}, (32) 

by employing the expression in Eq. (23), from Eq. (30) 
we obtain that (as in Eq. (23), below x ja): 



N/2 

G(fc; z) = a ^ e-*'^'(^") | e*'^'(^°) Gi,^i,(x, z) + e*'^'(^°) G^^i{x, z) + e^^'^'^^") G^^^ix, z) + e^''^^^^''^) G£,£(x, z)} 

i=-iV/2+l 
r.L/2 

Ax e-''^^ 

'-L/2 



"{e'^'^^ Gft.g(x,z)+e*- Gs.£(x, z) + e"*- G£,^(a;, ^) + e"*- G£,£(x,z)} 
G'R,R(fc - fcp; z) + GR_L(fc - fcp; z) + GL,R(fc + z) + GLx(fc + ^f; z) 
-|GR,R(fc - ^f; - z) + GR.L(fc - fcp; - ^) + Gi ^{k + fcp; - z) + GL.L(fc + ^f; - z)|, 

u - ■ u u - ■ U J 

I 



(33) 



where L and R are the time-scaled counterparts of L and 
R respectively, Eq. (C234), and u is the renormalized 
Fermi velocity Wp, Eq. (C123); in the limit V = Vc, one 
has u = ttvf/2 « 1.57 Eqs. (C40) and (C126). In 
arriving at the last equality in Eq. (33) , we have made use 
of the equality in Eq. (C235). The asymptotic equality 
'~' in Eq. (33) applies in the limit of a ^ 0, where, with 
L = Na, according to the Euler-Maclaurin summation 
formula one has (item 23.1.30 in Ref. [23]) 



N/2 
] = -N/2+l 



1 



rL/2 

a J-L/2+a 
1 



a 



L/2 



dx J{x) 
dx f{x). 



(34) 



In appendix C we deal with Gx3(fc;z), where A and B 

stand for one of L and R. 

In Fig. 2 we present the Green fimction G(fc; uo + iO+) 
for fc = fcp as calculated numerically on the basis of the 
expression Eq. (33), in which {Gx^^ik^k^; z/w)} are eval- 
uated according to Eqs. (C239) and (C272). We point out 
that the nature of the approximations underlying our cal- 
culations (specifically, the approximation of Eq. (C249)) 
restrict the validity of the calculated G(fc; w -|- iO"*") to the 
regions of the k-uj plane close to (±fcF,/i), where ^ = 
in our calculations. The existence of gap in the single- 
particle spectrum of the system under consideration, as 
evidenced by Im[G(fcF;'i^ + iO+)] = for |w|/(mu) < 2, 



trivially leads to the satisfaction of the Luttinger theo- 
rem, as discussed in Sec. IV. With reference to Eq. (22), 
for the mass gap one has 



7W = 2um, 



(35) 



where m is the soliton/antisoliton mass introduced in 
Eq. (C248); from Eq. (35) and Eq. (22), one readily de- 
duces the expression for m specific to the case at hand. 
Fig. 2 further reveals that our calculated Green function 
exactly satisfies Eq. (9) for fc = fcp and = 0, whereby 
fct is to be identified with fcp = 7r/2. This result is a di- 
rect consequence of the symmetry relation in Eq. (C241), 
which, on the condition of continuity of G{kp; z) at z = 0, 
leads to Eq. (C242). 



For completeness, the contributions to G(fc'F 



as depicted in Fig. 2, are essentially due to GR_R(0;a; -I- 
iO+) and Gr_l(0;w -|- iO+), Eq. (33), which are moreover 
equal. 



B. The case of F > t « t' 

In the limit of infinitely large V, one of the two sites 
in the unit cell of the CDW GS is fully occupied and the 
other is fully empty. Thus ioi V t t' , the effec- 
tive tunneling matrix element corresponding to nearest- 
neighbours is proportional to t'^/V, for which one has 
t'^/V ^ t' . It follows that in the case at hand, one can 
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FIG. 2: (Colour online) The real and imaginary parts of 
G(A;; + 10"*") for k — = 'k/2 as function of w, calculated on 
the basis of the expression in Eq. (33), in which {Gf,^s{k; iu^n)} 
are evaluated according to Eqs. (C239) and (C272); follow- 
ing Eq. (C241), G{kT,;ijo — iG^) is for < exactly equal to 
— G{kr; —ui + iO"*"). We have employed £, — 50, which, follow- 
ing Eq. (C253), corresponds to K = 25/51 < 1/2, for which 
the GS is insulating (see the remarks centred at Eq. (C97), 
the second paragraph of Sec. C4 and Sec. C4c). We have 
chosen a = 1, t — 1 and arbitrarily equated m with f/1000. 
Further, we have identified u with the value corresponding 
to A" = 1/2, namely Eq. (C126), where v-e — 2at, 

Eq. (C40). The overall magnitude of lm[G{k;Lo + iO+)] di- 
minishes rapidly for k departing from fcp (this is accompanied 
by an increasingly sharper rise of the right-most flank of the 
deep in Ke[G{k;uj + iO'^)] and a shift towards lower energies 
of the zero in this function). This signifies that the farther 
removed k is from ftp , the more significant become the higher- 
order soliton/anti-soliton contributions to G{k\LO + iO"'') that 
are neglected in the approximation of Eq. (C249). As ex- 
pected (since the minimum value of x, Eq. (C279), amounts 
to 2m), the threshold at which Im[G(fcF;a; -I- iO'*')] ^ satis- 
fies |Lj|/(mu) = 2. Note that G(A:f;0) = 0, which, in view of 
Eq. (C242), is an exact result. This outcome is significant in 
that in the present case A^l = ftp = 7r/2 and ^ = Q (compare 
with Eq. (9)). 

to a good approximation identify t with zero and only 
consider next-nearest-neighbour hopping processes. As- 
suming that for V = oo the odd sites are fully occupied, 
and thus the even sites are fully empty, for y ^ i w 
to a good approximation one has that electrons hop be- 
tween even sites and holes hop between odd sites. One 
thus trivially obtains from Eq. (4) that in the case at 
hand the energy dispersions of the single-particle exci- 
tations of the system under consideration are to leading 
order described by 

i?±~-2t'cos(2fc)+| ^ , fcG(-|,|]. (36) 



In the present case, the minimum spectral gap is an in- 
direct one and its value amounts to ^ — At' . 

VI. CONCLUDING REMARKS 

We have unequivocally demonstrated that the finite- 
size-scaling analysis as employed by KP [12], which is 
based on the explicitly-calculated values of the G{k\ ji) 
pertaining to finite systems consisting of N lattice sites, 
with N taking the values presented in Eq. (8), is not 
capable of establishing breakdown of the Luttinger theo- 
rem for the insulating uniform GSs of the t-t'-V model; 
using a simple model self-energy for which the corre- 
sponding G(fc; /i) exactly satisfies the Luttinger theorem, 
the approach by KP [12] implies failure of the Luttinger 
model. Furthermore, we have shown that the statement 
by KP [12], that the numerical results corresponding to 
A'^ = 26 implied breakdown of the Luttinger theorem, is 
erroneous; through a simple counting of the number of k 
points for which G{k; n) > (see Fig. 2 in Ref. [12]), one 
readily ascertains that for A^ = 26, which according to 
Eq. (7) corresponds to A'o = 13 at half-filling, the Lut- 
tinger theorem indeed applies: for exactly 13 different k 
points, out of the total of 26 different k points of which 
the underlying IBZ consists, Eq. (6), one has G{k; fi) > 0. 

On the basis of their calculations corresponding to the 
insulating GSs of the Hubbard Hamiltonian on a two- 
dimensional triangular lattice, KP in an earlier paper 
[24] arrived at a similar conclusion as the above con- 
cerning the Luttinger theorem, namely that in the ab- 
sence of particle-hole symmetry the Luttinger theorem 
fails. In a separate Comment [25], one of us unequivo- 
cally demonstrates that this conclusion is also false and 
that the quantitative amount of the reported "failure" 
of the Luttinger theorem is fully attributable to an er- 
ror arising from truncating by KP to leading order in 
p, = n — U/2 the infinite series expansion for G{k; n) in 
powers of jl that has underlain the considerations by KP 
[24]. 
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APPENDIX A: SOME COMPUTATIONAL DETAILS 



We have calculated the numerical results presented in Table I with the aid of Mathematica (version 6). The reader 
wishing to carry out calculations for the values of the parameters {/x, t, A, B} different from those to which the data 
in Table I correspond, may find the following Mathematica instructions helpful. The numerical results obtained for 
the specific values of fi, t, A and B that we choose below, are presented in the fourth row of Table I. As in Ref. [12], 
for N we have chosen the five values listed in Eq. (8). 

In[l]:= Ek[k_, t_, A_, B_] := -2 t Cos [k] + A + B/(.k'2 - (Pi/2)-2) 
In[2]:= G[k_, niu_, t_, A_, B_] := l/(mu - Ek[k, t. A, B] ) 

In[3]:= Gbar[mu_, t_, A_, B_, n_] := (G [Pi/2 + Pi/n, mu, t, A, B] + G[Pl/2 - Pl/n, mu,t,A, B])/2 

In[4]:= DeltaG[mu_, t_, A_, B_, n_] := G[Pi/2 + Pi/n, mu, t. A, B] - G[Pi/2 - Pi/n, mu, t. A, B] 

In[5]:= TGbar = Table [{n,Gbar [0 . , 1.0, 1.87, 1.33, n]}, {n, {14, 18, 22, 26, 30}}] 

Out [5]:= {{14, 0.701368}, {18, 0.367454}, {22, 0.228655}, {26, 0.156993}, {30, 0.114867}} 

In[6]:= InputForm [FindFit [TGbar , al + bl/n + cl/n"2, {al , bl, cl}, n] , NumberMarks -> False] 

Out [6] //InputForm= {al -> 0.19006514037109792, bl -> -10.33960282911629, cl -> 244.78744001923275} 

In[7]:= TDeltaG = Table[{n, DeltaG[0.0, 1.0, 1.87, 1.33, n]}, {n, {14, 18, 22, 26, 30}}] 

Out[7]:= {{14, -1.89343}, {18, -1.17785}, {22, -0.85808}, {26, -0.678616}, {30, -0.563547}} 

In[8]:= InputForm [FindFit [TDeltaG, b2/n + c2/n"2, {b2, c2}, n] , NumberMarks -> False] 

Out [8] //InputForm= {b2 -> -6.989008699335287, c2 -> -269.1504632514183} 

In [9]:= N[l/2 - 2 (0 . 19006514037109792/-6 . 989008699335287) , 7] 

Out [9] := 0.55439 



APPENDIX B: THE SINGLE-PARTICLE GREEN 
MATRIX OF THE CDW STATE 

Here wc determine the general form of the single- 
particle Green matrix G(fc; z) pertaining to the CDW 
state of the one-dimensional t-t'-V model that we con- 
sider in Sec. IV. Amongst others, we show that this 
matrix is a 2 x 2 one. As we shall see, this result is a 
direct consequence of the fact that the model under con- 
sideration is defined on a lattice; had the model been a 
continuum model, consisting of N/2 unit cells (each con- 
taining two 'atoms'), the relevant G{k,z) would be an 
N/2 X N/2 matrix [26]. 

For the regular lattice under consideration, with lat- 
tice constant a, let {\xj) \\j = l,2,...,iV} denote the 
complete set of the normalized eigenkcts of the position 
operators x, and {|fc) \\k E IBZ} the complete set of 
the normalized eigenkets of the momentum operator p. 
For the following considerations it is more convenient to 
renumber the lattice sites in such a way that 

r N N N^ 

je{-Y + l,-y + 2,... ,0,1,. ..,-}, (Bl) 

where N is even. Note that since the system under con- 
sideration is in the thermodynamic limit. N is macro- 
scopically large. One has 

(a;,|fc) = ^ e'*^^^''^) . (B2) 

^ Vn 

In the following, and '^j, denote sums in which j 
and j' traverse the set in Eq. (Bl). 

With G(fc, k'; z) denoting the momentum representa- 
tion of the Green operator G(z), of which Gjj'{z) is the 
site representation, by the completeness of the set 



one has 

G{k,k';z) =^{k\xj)Gj^f{z){xj>\k'). (B3) 
id' 

On account of the fact that in the CDW state the unit 
cell consists of 2 lattice points, one has the following sym- 
metry relationship: 

Gj+2^,f+2Uz) = Gjj'(z), (B4) 

for any finite integer value of v {i> may be negative, zero 
and positive). Using Eq. (B4), from the expression in 
Eq. (B3) one trivially obtains that 

G(fc, k';z) = c-2''i(''^-'^')° G(fc, k'; z), (B5) 

from which one deduces that unless 

k-k'^-e = K'f, ^ = 0,±1,..., (B6) 
a 

one has G{k,k';z) = 0. One observes that {K'^} is the 
set of the reciprocal-lattice vectors corresponding to the 
CDW GS of the system under consideration. Since k,k' G 
IBZ, it follows that the £ in Eq. (B6) can take only the 
three values —1,0 and +1. With 

1BZ' = (- — ,— ], (B7) 
^ 2a' 2a^' ^ ' 

below we shall denote all non-trivial Fourier components 
of G(fc, k'; z) by Gi^c {k; z), where 

Gi^i>{k;z)^G{k + K'i,k + K'e,;z), fc e IBZ'. (B8) 

Following Eq. (B3), one thus has (see Eqs. (B2) and 
(B6)): 

Ge,Ak; z) = Y.^k + K[\x,)G,.y{z){x,, \k + K',,) 

3,3' 

^ ^(-1)^-^+/^' {k\x,)G,,,,{z){xr\k). (B9) 

3,3' 
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Using the closure relation for {\k) \\k E IBZ}, one can and further 
write 



(BIO) 



fc,fc'eiBZ 
Decomposing the IBZ as 

lBZ = (--,-f]UlBZ'u(f ,-], (Bll) 
a Za Za a 



IBZ' = (--^, 0] U (0, -^1 = 1BZ'_ U IBZ' (B12) 
2a 2a 



on account of the fact that the only non-trivial Fourier 
coefficients G{k,k';z) are those presented in Eq. (B8), 
after some algebra one obtains that 



J 



Gjj'iz) = {i^j\'^)'^o,o{k;z){k\xf) + (xj\k)Go,i{k; z){k + K[\xj>) + {xj\k + K[)Gi^o{k; z){k\xj>) 

fceiBZ'_ 

+ (xj\k + K[)Gi^i{k; z){k + K[\xj>)^ 
+ E [{xj\k)Go,o{k;z){k\xj,) + {x,\k)Go^^i{k;z){k-K[\xj,) + {x,\k-K[)G-i^^^^^^ 

+ (x^\k - K[)G^i,^i{k; z){k - K[\x,>)y 



fceiBZ 



(B13) 



In arriving at this expression we have used the fact that 
e.g. the interval (^f is translatable into IBZ^ by 

K[ . We note that the total number of terms contributing 
to Gjji (z) is nine, of which three are identically vanishing 
(such as the one corresponding to J2keiBZ' Sfc'eiez' ) 
on account of the fact that the pertinent k and k' do not 
satisfy the relationship in Eq. (B6). 
From Eq. (B2) one obtains that 



c,\k±K{) = {-l) 



i\k), 



(B14) 



where the RHS is independent of whether on the left- 
hand side (LHS) one has +K[ or —K[. Using the result 
in Eq. (B14), from Eq. (B3) one immediately deduces 
that 

G-i^-i>{k;z) = Ge,e'ik;z), £,£' e {-1,0,1}. (B15) 

Using Eqs. (B14) and (B15), Eq. (B13) can thus be ex- 
pressed as 



Gj,A^)= E v]{k)-G{k;z)-v,,{k), 
feeiBZ 



where 

Go,o{k;z) Go,i(fc;z) 
Gi.o{k;z) Gi^i{k;z) 
and (cf. Eqs. (B2) and (B14)) 



G(fc;z) 



(B16) 



(B17) 



vjik) 



{k + K[\x, 



N 




The vector Vj{k) is a row vector whose first and second 
components are {xj\k) and {xj\k + K[) respectively. 

That G(fc; z) is a finite-dimensional matrix, is related 
to the fact that the IBZ is finite, on account of the system 
under consideration being defined on a regular lattice of 
non-vanishing lattice constant; that it is a 2 x 2 matrix, 
is a consequence of the fact that the extent of IBZ is 
twice as large as that of IBZ'. Note that for the homo- 
geneous GS of the system under consideration (that is 
for V < Vc), Go,i(fc;z) = Gi,o(fc;^) = 0, Vfc e IBZ'. 
In such case, Go^o{k;z) and Gi^i(fc;z), with k G IBZ', 
together describe the Green function G{k; z) of the uni- 
form GS, with k S IBZ. In this connection, on account of 
Eq. (B15), Gi i(fc; z) for k e IBZ'^ is identical to G{k; z) 

forfce(-f,-i,]. 



APPENDIX C: THE BOSONIZATION OF THE 
t-t'-V MODEL, THE QUANTUM SINE-GORDON 
AND THE XXZ SPIN-CHAIN MODELS 

Here we determine the bosonized version of the one- 
dimensional t-t'-V model in the continuum limit [27], 
where (below a the lattice constant and j a site index) 



- JO. 



(CI) 



is considered as a continuous variable. In doing so, we 
formally assume that a ~* and A'^ — > oo, in such a way 
that 



L = Na 



(C2) 



(B18) 



is macroscopically large. In practice, the lattice constant 
a is identified as the unit of length so that physically a — s- 
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is to be understood as signifying that the correlation 
functions calculated in the continuum limit accurately 
describe those pertaining to the underlying lattice model 
over the length scales that are large in comparison with 
a. 

Unless we indicate otherwise, in the following we shall 
explicitly deal with half-filled GSs. for which the Fermi 
points of the underlying non-interacting GS is located at 

±k, = ±^. (C3) 
Za 



1. Preliminaries 

Making use of the exact representation of Cj in terms 
of its Fourier components {ofe}, that is (cf. Eq. (B2)) 

feeiBZ 

disregarding the contributions corresponding to values of 
k outside the interval 



[-fcp - A, -fcp + A]\J l+kp ~ A, +kF + A]<Z IBZ, 

where yl > 0, one obtains that 
c,c. V^{c+'^^(^''' ^ ^ e'^^f^-^) afc+,,+e-* ^ e'"^'"^ dk-k^^V^{{+^yRix) + {-^yL{x)}, (C5) 

\k\<A |^,|<^ 

where the (iz)-' are the values of e**'^'''^-'''^ for the specific k^ in Eq. (C3). 

The slowly-varying chiral field operators R{x) and L{x) can be expressed in terms of the bosonic field operators 
S{x) and 0{x) as follows (cf. Eq. (2.30) in Ref. [28]) [29]: 

1 X e^*'^-" = lim ^±^{k.-^/L)x exp[±zV^<^(a;) - iV^0{x)] 

= e=F*''/^ lim g±»(fcp-,r/L)x g^pF i^0[x)] exp[±iV7r<^(x + 0+)] , (C6) 

"^0 V27ra 



where tj± is the Klein factor [19, 28] and a « 1/yl. Nat- 
urally, the operators (I>{x) and 0{x) arc functions of this 
cut-off parameter. Since the problem at hand is defined 
on a lattice, of lattice constant a, it is natural to identify 
a with a. This is the choice that we shall make in the 
following. 

The second equality in Eq. (C6) is established through 
the application of the Baker-Campbell-Hausdorff formula 
[30] 

gAgB^gi+Bgi[A,Bl-^ (C7) 

which applies when [A, B]_ commutes with both A and 
B, and making use of Eq. (CIS) below. With reference 
to the latter equation, it is evident that the phase factors 
g=F!7r/4 ^Yic second expression in Eq. (C6) are directly 



related to our use of 'P{x + 0+) {and L oo), instead 
of •P{x). Use of <P{x + 0+) is fundamental, as it pre- 
vents the possibility that dx0{x) may occur alongside 
(l>{x) (or 0{x) alongside dx'i'{x)) whose commutation is 
unbounded; see Eq. (24) above. We note in passing that 
the first expression in Eq. (C6) coincides with those in 
Rcfs. [31], [4] and [28]. As we shall see later (Sec. C6), 
the second expression in Eq. (C6) will he crucial in es- 
tablishing a relationship between R{x) and L{x) and the 
soliton-generating fields pertaining to the quantum sine- 
Gordon Hamiltonian [21] whose form factors have been 
determined by LZ [32]. 

For d^(x) and 0{x) one has the explicit expressions (cf. 
Eq. (2.23) in Ref. [28]; [29]) 



${x) 



^X ^ , «V7F ^i|fc|^l/2 1 ^-a\k\l2-tkxr^\ 

V^.T - Ja/tF /^Ifcls 1/2 1 



L 2^ ' 



k' 



0ix) = -^i-^E(^) '^^-""^'/^-"^(^I-M, 



(C8) 
(C9) 
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in which and {bk} are canonical bosonic creation 
and annihilation operators, respectively. Above 

Q = Qr + Ol (CIO) 

is the total-charge operator, and 

J = Qn-Qi. (Cll) 

the total-current operator. Since Q and J commute with 
the Hamiltonian under consideration, the Hilbert space 
of the problem at hand can be classified in terms of the 
eigenvalues Q and J of these operators, which are integer 
valued and satisfy the condition [19, 33] (see also p. 280 
in Ref. [20]) 

(-1)-^ = (-l)Q. (C12) 

For clarity, our Q is in Ref. [19] (in contrast to Refs. [33] 
and [34]) denoted by iV — Nq, where A'^o = fepi/Tr is as- 
sumed to be an odd integer. This fact accounts for the 
absence of a minus sign on the RHS of Eq. (C12) in com- 
parison with Eq. (2.3) in Ref. [19]. 

For the explicit calculations to be carried out below, 
it is relevant to note that on the basis of the expressions 
in Eq. (C8) and (C9) one can easily demonstrate that 
(below X, x' are arbitrary) 

(i) ^[x) and 0{x) are Hermitian, 

(ii) «?(a;) commutes with ${x'), and 

(iii) 0{x) commutes with 0{x'). 

Since x and x' are arbitrary, from (ii) and (iii) one im- 
mediately infers that 

(iv) dx'P{x) commutes with ^{x'), and 

(v) dxO{x) commutes with 0{x'). 

On account of (ii) and (iv) , in dealing with the functions 
of <P{x) one can proceed as though 'l>{x) were an ordinary 
function, and not an operator. Similarly for 0{x), on 
account of (iii) and (v). In this connection, and in view 
of the exponential functions of (!>{x) and 0{x) that we 
shall encounter below, we point out that for an arbitrary 
operator A{x) one has [35] 

dx e^("' = ds [dxA{x)] e^^(") , (C13) 

which, for A{x) and dxA{x) not commuting, amounts to 
a non-trivial and generally a very complicated expression 
to work with. 

It is important to specify the nature of the sums with 
respect to k in Eqs. (C8) and (C9). These sums are over 
[19] 

|fc= ^||n = 0,±l,±2,...,±oo|. (C14) 



This set is different from the discrete set 

r . 27rn „ L L ^ s 

{^ = — "|-2^<"^2^} (^^^^ 

comprising the IBZ corresponding to the system under 
consideration, Eq. (6). Thus, insofar as the sums in 
Eqs. (C8) and (C9) (and those in the related expressions) 
are concerned, the usual sharp cut-off of summations of 
the fc-space, as imposed by the hounded domain IBZ, is 
replaced by the soft cut-off e""''''. This replacement is 
crucial for the validity of some of the expressions that 
we shall encounter below. 

From Eqs. (C8) and (C9) one readily deduces that 

^{x + L) = S{x) + ^Q, (C16) 
0ix + L) = 0{x)-^J. (C17) 

For completeness, the Qr and Ql in Q = Qr + Ql are 
in Ref. [28] denoted by N-^ and Ni^ and the eigenvalues 
corresponding to these by respectively M-n and M-l (see 
Appendix B in Ref. [28] and in particular note that Hal- 
dane's above-mentioned TVq is in this Appendix denoted 
by 2no + 1). 

Making use of the explicit expressions in Eqs. (C8) and 
(C9), and the canonical commutations relations for the 
bosonic operators {b\^ and {6fe}, one trivially obtains 
that for a;,a;' G R (see Eq. (2.24) in Ref. [28]) 

[^{x), 0{x')\ _^]^^\ ^^k[x'-x)-^\k\ 

_i _i / c-2'^"/-^sin(27r[a;'-x]/L) \ 
TT VI- e-^'^a/i cos(27r[x' - x]/L) J 

1 

^ - sgn(x' — a;) when a I 0, L ^ oo. (C18) 

The last expression in Eq. (C18) is the leading-order 
asymptotic contribution (corresponding to the asymp- 
totic limit L -~^ oo) to the function on the second line, 
subject to the condition that \x' — x\ L; the expres- 
sion can be treated as an exact one (i.e. the '~' can 
be replaced by '=') by taking the limit L oo prior to 
allowing x' — a; to vary over the interval (— oo, cx)) . We re- 
mark that on employing the set of k points in Eq. (C15), 
instead of the one in Eq. (C14), one would obtain a func- 
tion that for ]a:' — x] ^ L is up to a finite multiplicative 
constant almost identical to Si(7r[a;' — x]), where Si(z) is 
the sine-integral function (item 5.2.1 in Ref. [23]); im- 
portantly, the derivative of this function with respect to 
x' is finite at x' = x for all values of L (compare with 
Eq. (24) in the light of Eq. (C19)). 

We emphasise that the commutator in Eq. (C18) is a 
periodic function of x'—x, with period L. Since for x, x' G 
[0, L] one has x' ~ x E [—L, L], it follows that unless one 
takes the limit L ^ oo prior to allowing x' — x to vary 
over the interval (—00,00), use of the last expression in 
Eq. (C18), which is not a periodic function of x'— x, is not 
justified. This fact has not be taken account of in Refs. [4, 
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28, 31]. In contrast, Haldane [19] takes careful account of 
the periodicity of the commutator in Eq. (CIS), as well 
as other related commutators. In doing so. Haldane only 
considers the limit a = 0"*", corresponding to the case 
where the cut-ofF function e~"l'"'l is infinitely soft. 
Taking the limit L oo, for 

nix) = d.eix) (C19) 

from Eq. (C18) one obtains the commutation relation in 
Eq. (24). It follows that for L = oo, n{x) is the mo- 
mentum operator canonically conjugate to 0{x). With 
reference to our earlier remarks, for finite but sufficiently 
large values of L, instead of Eq. (24) one has the following 
asymptotic expression (cf. Eq. (3.27) in Ref. [19]) 

oo 

[S{x),n{x')]_^i d{x-x'+nL), 

n— — oo 

(C20) 

Within the framework of the bosonization technique, 
one defines the chiral density-fluctuation operators 

Mx) = :R\x)R{x):, p^{x) ^ -.L^ {x)L{x) (C21) 

which are related to the bosonic fields 'P{x) and 0{x) as 
follows (cf. Eq. (C19)): 

dJix) = {Mx) + Pl{x)) = V^p{x), (C22) 

d,0{x) = -V^ip^x)^ p^x)) =n{x). (C23) 

In Eq. (C21). -.AB: stands for the normal-ordered coun- 
terpart of the product AB. The expressions in Eqs. (C22) 
and (C23) are readily obtained from Eqs. (C8) and (C9) 
by employing the expression for the canonical boson op- 
erator bk in terms of the Fourier components of the chi- 
ral density operators Pr{x) and pi^{x) (sec Eq. (2.15) in 
Ref. [28]). 

We note that pjt{x) + pi^(x) is the slowly-varying part of 
the total density-fluctuation operator [36] . To appreciate 
this fact, one has to realise that, making use of Eq. (C5), 

« ■.R\ja)RUa): + : L\ja)Lija): 

a 

+ {-iy [■.R\ja)L{ja): + : {ja)R{ja) :]. (C24) 

One observes that Pr{x) + pi^{x) does not take account of 
the last two terms on the RHS of Eq. (C24). As we shall 
see later (see e.g. Eq. (C31) below), the last-mentioned 
two terms do not contribute to the total number of parti- 
cles in the system to leading order in the small parameter 
a. From Eqs. (C16) - (C22) one has 

fL ( Q, 

J dx [p4x) ± Pl(x)] I ^ (C25) 

These expressions make explicit that Q and J are the k = 
Fourier components of Pr{x) + Pi^{x) and Pr{x) — Phi^) 
respectively. 



To make contact with the work by LZ [32], we now 
proceed with expressing 0{x) in terms of IJ{x). To do 
so, it is required that 0{x) be specified at one value of x; 
following Eq. (C19), II{x) determines 0{x) up to an ad- 
ditive constant. With reference to the Riemann-Lebesgue 
lemma (Sec. 9.41 in Ref. [37]), the second terms on the 
RHS of Eq. (C9) becomes vanishingly small for |a;| oo, 
whereby 

0{-L)^^J for L^oo. (C26) 

Thus, integrating both sides of Eq. (C19) with respect to 
X, one obtains that [38] 

0{x)^^/^J + j dx' n{x') for L^oo, (C27) 

where '~' is to be identified with '=' for L = oo. In a 
Hilbert space corresponding to a fixed eigenvalue J of J, 
the operator y/nJ on the RHS of Eq. (C27) gives rise to 
a constant phase factor in the expressions for R{x) and 
L{x), which, with reference to Eq. (C6) and in view of the 
fact that J takes integer values, is equal to (— 1)"^ ~ ±1, 
depending on whether J is even or odd. This factor is of 
no consequence to $^a.b(3^i ''"); ^^^d the related functions, 
when A ^ B, where A,Be {L, R}; it is however relevant 
when A^ 13. To clarify this, let \W) denote an arbitrary 
eigenstate of Qr and Ql in the Fock space of the system 
under consideration. With reference to Eq. (Cll), and 
since the J in J|'f') = is the difference in the num- 
bers of right-moving and left-moving particles in \^), in 
addition to the fact that the number of the left-moving 
(right-moving) particles in i^ltf') (in W\9)) is by one unit 
more than that in \^), one immediately infers that the 
expectation value of RU (or of Lw) with respect to \^) 
differs by a minus sign from the one corresponding to the 
case in which the -^/tF J on the RHS of Eq. (C27) has been 
suppressed altogether. The '±' pre- multiplying the ex- 
pression on the RHS of Eq. (25) (which we could equally 
have chosen to be '=f') accounts for the last-mentioned 
minus sign. As we shall see later in this appendix, with- 
out this minus sign, the contributions of G^ i^k — k^; z) 
and Gi^ii{k + k^iz) to the calculated G{k;z), Eq. (33), 
would violate the causality principle (see in particular 
the second remark following Eq. (C236) below). 

For L — *■ oo, one generally identifies the Q^^i-^/^)^ on 
the RHS of Eq. (C6) with unity. Thus by equating a 
with the lattice constant a (see above) and making use 
of Eq. (C27), for L ^ oo the expressions in Eq. (C6) 
are transformed into the those presented in Eq. (25). In 
our calculations we employ the expressions in Eq. (25) by 
further suppressing the Klein factors U± , this on account 
of the fact that U± do not affect the behaviours of corre- 
lation functions that we consider in this paper [39] (for a 
case where Klein factors play a role, consult Sec. 4.3.2 in 
Ref. [28]). Consequently, the expressions in Eq. (25) in 
which the U± are suppressed, are not operator identities. 
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however can be viewed as such in calculating the cor- 
relation functions {Gf^^^{x,iujn)} introduced in Eq. (31). 
Barring the functions e^''^^^-'"' that we have already sep- 
arated from R{x) and L{x), the expression on the RHS 
of Eq. (25) concerning _R(x) {L{x)) coincides with that of 
the operator Oi (O2) introduced by Luther and Pcschel 
in Ref. [39]. 

We note in passing that, with II {x) being associated 
with the quantum- mechanical momentum operator p, the 
contribution to R{x) and L(x) of the first term in the 
argument of the exponential function on the RHS of 
Eq. (25) is the equivalent of the single-particle eigenfunc- 
tion exp{ipx) of the one-dimensional Schrodinger equa- 
tion describing a particle in a constant potential. Fur- 
ther, from Eq. (C22) one observes that placing a particle 
at location x = xo gives rise to a 'kink' (a discontinu- 
ity if the particle is infinitely localised), of magnitude 
y/n, in <P{x) at a; = xq. Therefore, physically the fac- 



J 



tor ex.p[±i^ S{x)] on the RHS of Eq. (25) introduces a 
phase jump equal to -I-tt (— tt) in R(x) {L{x)) for any x 
at which a localised particle is encountered [19, 40] (see 
in particular Sec. 3.1, "Phenomenological bosonization" , 
in Rcf. [28]). The accumulation of these phase jumps 
amount to the Jordan- Wigner phase that ensures anti- 
commutation relations to be satisfied for Fermi operators 
expressed in terms of Bose operators. 

The dynamics of the bosonic field operators <P{x) 
and II(x) are determined by the quantum sine-Gordon 
Hamiltonian [21]. Below we present the details un- 
derlying the derivation of this Hamiltonian for the 
one-dimensional t-t'-V Hamiltonian under consideration. 
This derivation almost entirely coincides with that en- 
countered in the treatments of the Heisenberg-Ising 
chain of the SU(2) spin operators, as described in e.g. 
Refs. [4, 28, 31, 41]. 



Bosonization and tlie sine-Gordon Hamiltonian 



Using the expression in Eq. (C5), one trivially obtains that 

c]cj+i ~ ai{RHx)R{x + a) - D{x)L{x + a) - {-ly R'^{x)L{x + a) + {-ly L'^{x)R{x + a)}. (C28) 
Making use of the Eulcr-Maclaurin summation formula (item 23.1.30 in Rcf. [23]), whereby 



^fUa)^- dx 



/(^) 



(C29) 



for a — !■ one obtains that 

N 



RHja)R{ja + a) c^- [ dx R\x) R{x + a) ^ - [ dx R\x)R{x) + [ dx R'' {x)d^R{x). 
~{ a Jo a Jo Jo 



(C30) 



For X]j=i ^'\i^)L{jo. + a) one obtains a similar expression. 

For obtaining the sum with respect to j of the last two terms on the RHS of Eq. (C28), we assume that N is even, 
whereby one can write 

AT N/2 N/2 

Y^{~iy R^{x)L{x + a) = ^i?tQ-(2a))L(j(2a) + a) - ^i?t(j(2a) - a)L{j{2a)) 
j=i j=i j=i 

wi/ dx R^{x)d^L{x) + ]-( dx{d^R\x))L{x):^l{Rl{L)L{L)^R^{2a)L{2a))^0{a), (C31) 



where the last equality follows from the periodic bound- 
ary condition 

R'<{x + L)L{x + L) = R\x)L{x). (C32) 

This condition can be readily deduced from Eq. (C68) 
below in conjunction with Eq. (C16) above, taking into 
account that the eigenvalues of Q are integers. 



For J2jLii^^y + oiic obtains a similar 

result, implying that for a ^ the last two terms on 
the RHS of Eq. (C28) do not contribute to 
This is a direct consequence of half-filling for which kp 
has the value given in Eq. (C3). 

Combining the above results, for a ^ and N ^ oo 
(see Eq. (C2)) one obtains that 
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In a similar fashion as above, one obtains that 

N 



i( Vc]cj+i +H.C.) = -{2at) i / dx {R\x)d^R{x) - D {x)d.^L{x)) . 



Of' 

t'( Vc]cj+2 +H.C.) = — / dx {R^x)R{x)+D{x)Lix)). 
3=1 ^ 



(C33) 



(C34) 



In view of Eqs. (C21), (C22) and (C16) one has 



at; 



N 

^:c;.c,+2: + 



H.c.) = 



2t' 

-Q 

a 



dx dx^{x) 



(C35) 



This term can be absorbed in the contribution of the 
chemical potential in the thermodynamic Hamiltonian 
of the problem at hand (see in particular Eq. (C25) 
above). It follows that at half-filling^ and for sufficiently 
small V/t, whereby the bosonization scheme as consid- 
ered above is physically reasonable, the parameter t' does 
not enter into the bosonized Hamiltonian. As we shall 
discuss later, this Hamiltonian only to leading order in 
V/t coincides with the bosonized Hamiltonian derived 
from the exact solution of the t-V Hamiltonian (compare 
the n in Eq. (C75) with that in Eq. (C76)). 

For completeness, we point out that the results in 
Eqs. (C33) and (C35) are immediately obtained by em- 
ploying the following exact representation (see Eq. (C4)): 

N 

-^{tc]cj + i+t' c]cj+2)^^-C-^ X] 

j=i fceiBZ 

(C36) 

where for the energy dispersion e{k) one has 

e{k) = -2t cos(fca) - 2t' cos{2ka). (C37) 
On effecting (cf. Eq. (C5)) 

^ e{k)alak^ ^ e(fc)a^afe, (C38) 



fceiBZ 



'A<k+k^<A 
'A<k-k^<A 



where ^ > is the cut-off in the momentum space, and 
linearising e(fc) at fc = ztfep, one has 



e{k) K -2{t + t')±v4kTTT/a), (C39) 



where 



Up = 2at, 



(C40) 



one obtains that 



e{k)alak -2{t + t') ^ a 



ak 



fceiBZ 



-A<k+kr<A 
-A<k-ki<A 



-Wf| X (fc-fcp)a^afc- ^ {k + kp) cildkY 

-A<k+kp<A 

(C41) 



-A<k-kp<A 



Evidently, Vp is independent of t' . The first term on the 
RHS of Eq. (C41) can be absorbed in the contribution of 
the chemical potential in the thermodynamic Hamilto- 
nian. With reference to Eq. (C40), it should be evident 
that the normal-ordered form of the expression on the 
RHS of Eq. (C33) coincides with the contribution pro- 
portional to Vf in Eq. (C41). 

We point out that the non-linear corrections to e(fc), 
with e{k) as expanded around k = ifcp, depend non- 
trivially on t' . However, since the contributions arising 
from these corrections correspond at the lowest order to 
four-fermion interactions in the Hamiltonian, they are 
all irrelevant operators, as can be established by simple 
power counting (for a precise treatment, one may consult 
Sec. 4.3 in Ref. [20]). It follows that for sufficiently small 
values of t', at half-filling the consequences oi t' ^ can 
be accounted for by perturbation theory. In our explicit 
calculations in this paper, we assume that t' is sufficiently 
small and thus neglect its consequences altogether. 

Following Eq. (C6) and the remarks regarding e.g. U± 
in the paragraph subsequent to Eq. (C27) on page 14, we 
write (cf. Eq. (14.31) in Ref. [20]) 



R{x) 

Ux) 



/27ra 



c-x^[-i^/^T {^^{x) + 0{x))\ 



\J2'Ka 



cxp[-iV7r6'(a;)] exp[±i\/7r ^(a; -I- 0+)], 



(C42) 



from which one deduces that (see the remarks (i)-(v) on 
page 13) 



dxR{x) 



-iMd.Oix) T dxi>{x)) { f J""^)' (C43) 

L{x). 



On account of the expressions in Eqs. (C21), (C22) and 
(C23), one thus obtains that (see Eqs. (C33), (C19) and 
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(C40)) 



(2at) i / da; {R\x)d^R{x) - L\x)d^L{x)) 



dx + {dj{x))^). (C44) 

Jo 

The right-hand side of this expression, which coincides 
with the expression for the Hamiltonian in Eq. (3.58) of 
Ref. [19], is twice as large as the correct non-interacting 
Luttinger Hamiltonian. The truth of this statement is 
directly verified by substituting the expressions for <P{x) 
and 0(2;) in Eqs. (C8) and (C9) in the expression on the 
RHS of Eq. (C44); following some algebra, one readily 
deduces that the resulting Hamiltonian is indeed twice 
the non-interacting Hamiltonian as expressed in terms of 
the bosonic operators {bk} and {b\}, such as presented 
in e.g. Eq. (2.22) of Ref. [28]. 

To circumvent the above-mentioned problem, which 
has not been recognised earlier, following Refs. [31] and 
[4] we proceed as follows. In order to keep the notation 
light, below we suppress the 0+ in ${x -\- 0"*"), with the 
understanding that when necessary the x in ^{x), and 
related operators, stands for x + . 

With 



it 



(C45) 



the imaginary time, we define the complex 'time' coordi- 
nates 



ixjvY 



ix/v^, 



(C46) 



and express the imaginary-time Heisenberg-picture coun- 
terparts of the operators ^{x) and 0{x), that is 

<^(a;,-ir) =<^(z,z), 6>(a;, -it) = 6>(z, z), (C47) 

in terms of respectively the analytic and anti-analytic 
operators (p{z) and (p{z), as follows (Eqs. (22.11) and 
(22.13) in Ref. [4]): 



0{z,-z) 



V'{z) + 'p[z), 



(C48) 
(C49) 



From the first expression in Eq. (C42), for the imaginary- 
time Heisenberg-picture counterparts of R{x) and L{x) 
one thus obtains that 



R{x, —ir) 
L{x, —ir) 



1 



'27ra 
1 

'27: a 



exp[+iV47r(^(z)], (C50) 



cxp[-iv47r(^(z)]. (C51) 



For R'^{x,-iT) {D{x,-iT)) one has the same expression 
as R{x,—iT) (L(x, —ir)) with -\-i {—i) however replaced 
by —i (+*)■ Compare with the relevant expressions in 
the bosonization tables 3.1 and 22.1 of Refs. [31] and [4], 
respectively. 



With f{z) and g{z) differcntiable functions of z and z, 
one has drf{z) = —ivpdxf{z) and drg{z) = -\-ivpdxg{z), 
so that 

dz^^idr-ivpOx), dz^^{d-r+ivpdx), (C52) 



and consequently 



dx = — {dz - dz) 



(C53) 



From Eqs. (C48) and (C49) one trivially obtains that 



dzSiz,z) = -d,iS{z,z) + diz,z)), (C54) 

dzS{z,z) = ^d-ASiz,z)^e{z,z)). (C55) 

One can thus make the following identifications (see the 
bosonization tables 3.1 and 22.1 of Refs. [31] and [4], 
respectively) : 



w{x, —iT)R{x, —ir) 
L\x, —iT)L(x, ^ir) 



+i 



-dz${z,z), (C56) 
-dzS{z,z). (C57) 



Making use of the equalities in Eqs. (C54) and (C55), 
one readily verifies that for t = the expressions 
in Eqs. (C56) and (C57) exactly reproduce Pr{x) and 
Pl{x) respectively, where Pr{x) and Pl{x) are defined in 
Eq. (C21) (see Eqs. (C22) and (C23)). 

In order to capture the correct prc-factor i in the 
bosonized non-interacting Hamiltonian, with reference to 
Eq. (C53) we propose the following substitutions in going 
from the LHS to the RHS of Eq. (C44): 

) ^ R''{x,0)dzR{x,0), 

D{x)dxL{x) ^ +—L\x,Q)dzL{x,Q). (C58) 

Vf 

These substitutions are motivated by the fact that 
R{x,—iT) depends only on z, Eq. (C50), and L(x,—iT) 
depends only on z, Eq. (C51). Using the expressions in 
Eq. (C52), one trivially obtains that 

~ —RHx,0)dzR{x,0) = \R\x)dxR(x), 
Wf 2 

+ -Lt(a;,0)9,L(.T,0) = \L\x)dxU^), (C59) 
Up 2 

where the dxR{x) and dxL{x) on the RHSs are to be 
expressed in terms of dxO{x) and dx^{x) on the basis 
of the expressions in Eq. (C43). The pre- factors ^ on 
the RHSs of the expressions in Eq. (C59) account for 
the aforementioned missing i on the RHS of Eq. (C44). 
These ^ pre-factors would turn into 1 on replacing the 
dzR{x,Q) and dzL{x,Q) in Eq. (C58) by respectively 
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dzR(x, —iT)\T=o and dzL{x, —iT)\r=o- Thus, the substi- 
tutions in Eq. (C58) prescribe a particular scheme, ac- 
cording to which the total variations of R{x, —ir) and 
L(x, —ir) at r = are to be identified with their partial 
variations along the x axis. 

We thus obtain that (cf. Eq. (C44)) 

-{2at) i I Ax {R^{x)d^R{x) - D{x)d:^L{x)) 
Jo 

^ {idMx))-' + idJix))-') 

EE I dx {n\x) + {d^S{x) f) , (C60) 

where in arriving at the last expression we have employed 
Eq. (C19). Consequently, from Eqs. (C33), (C35) and 
(C40) one obtains that 

N 

-Y.(*c]c, + l+t'clc,+2)+Ii.C. 
j=l 

^YJo ^^{n\^) + {dA^)y}, (C61) 

where denotes the association between the original 
Hamiltonian on the left and its bosonized counterpart on 
the right. 

Following Eq. (C5), for fij = c^Cj one has 

fij ~ a{R^{x)R{x) + D{x)L{x) 

+ i-iy [R^x)L{x) + L^ix)Rix)] }, (C62) 

which on normal-ordering can be expressed as 

hj c^a{p{x) + {-iy M{x)}, (C63) 

where p{x) is the total-density operator, Eq. (C22), and 

M{x) = ■.R^x)L{x): + : D {x)R{x) :, (C64) 

an operator that measures the amount of chiral symmetry 
breaking. With reference to Eq. (20), from Eq. (C63) one 

I 



observes that, at half- filling the GS expectation value of 
ap(x) per unit length, i.e. Q/N, is equal to 1/2, and 
that the GS expectation value of aM{x) per unit length 
is equal to 5n. In this connection, and with reference to 
Eqs. (C2), (C25) and (C29), we note that 

j = l "'O 

As we shall indicate in Sec. C4, 5n is the counterpart 
of the staggered magnetization in the Ising anti- 
ferromagnetic GS of the XX Z Heisenbcrg Hamiltonian 
for SU(2) spin operators. 

Making use of the expression in Eq. (C63), one obtains 
that 

hjfij+i ~ {p{x)p{x + a) — M{x)M{x 4- a) 

-{-ly [p{x)M{x + a) - M{x)p{x + a)] }. (C65) 

By the same reasoning as leading to the result in 
Eq. (C31), the contributions to of the two 

terms on the RHS pre-multiplied by (—1)-' can be ne- 
glected in the limit a ^ 0. In what follows we shall 
therefore only consider the contributions of p{x)p{x + a) 
and M{x)M{x + a) to the latter sum. 
Following Eq. (C22), for a ^ one has 

p{x)p{x -t- a) ~ p{x)p{x) = -{d,^x)f. (C66) 

TT 

In order to determine M{x)M{x -t- a), we make use 
of the expressions in Eqs. (C50) and (C51). In doing 
so, we employ the Baker-Campbell-Hausdorff formula in 
Eq. (C7) to express products of exponentials of ip{z) 
and (p{z) as an exponential of •p{z) -f (p{z) = <P(z,z), 
Eq. (C48). As we are dealing with r = 0, we need 
only to determine [tp{ix/vF),(p{—ix'/vF)]-- Since ^(z,^) 
{0{z, z)) coincides with S{x + 0^) {0{x)) for r = 0, from 
Eqs. (C48) and (C49) and the commutation relation in 
Eq. (C18) (i.e. the one corresponding to a J, 0, L — > oo) 
one readily obtains that 



[ipi-),0i )]_ = -['Pix + 0+) + Oix),<Pix' + 0+)~O{x')]_ ^ -{sgn{x~x' + 0+)~sgn{x~x' -0+)}. (C67) 



I 

Thus one has These expressions can be directly obtained by making 

use of the first equality in Eq. (C42). From Eqs. (C64) 
. and (C68) one deduces that 

R^x,0)L{x,0) ~ ^exp[-iVI^S{x)], 

ZTTa 

L\x,0)R{x,0) —^cxp[+iV^S{x)]. (C68) M(x) ~ — : sin(v^<P(x)) : . (C69) 

ZTTa Tra 
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We note that without the 0+ in ^{x + 0+) in the expres- 
sions for R{x) and L{x) (here, those according to the 
second equality in Eq. (C42)), the pre- factors in the ex- 
pressions in Eq. (C68) would have involved -1-1 instead of 
±i, resulting in M{x) ~ — — :cos{^/A^^^{x))^. (compare 
with Eqs. (6.23) and (6.24)Tn Ref. [28]; [29]). 

Since ^(x) commutes with 3{x + a) (see remark (ii) on 
page 13), in calculating M{x)M{x + a) one can consider 
^(x) and 3{x + a) as c numbers. Making use of the 
identity 



and of the expansion 



${x + a) c:i${x) +ad^${x) for a ^ 0, (C71) 



as well as the relationship 



:cos(i): « --i^, 



(C72) 



2 sin(.T) sin(y) = cos(a; — J/) — cos(.t -|- y), (C70) one readily obtains that 



Mix)M{x + a) 



-1 



2(7ra)2 

From Eqs. (C65), (C66) and (C73), one arrives at 

rL 



^S{x)) - ^{d^S{x))^. 



N 



f dx \ — {d^S{x)f + -^cos{VT6^S{x))} 



as a — > 0. 



(C73) 



(C74) 



Combining this result with that in Eq. (C61), one obtains the following correspondence between H and the bosonic 
Hamiltonian TL in the continuum limit: 



dx 



{|[iI^(x) + (l-K^)M(x))V 



2aV 
(2^ 



■ cos(v 167r 



S{x))}. 



(C75) 



Provided that V/t is sufficiently small, this sine-Gordon 
Hamiltonian [4, 21, 31] governs the dynamics of the field 
operators 11 (x) and ${x). With reference to our remark 
following Eq. (C69) above, we note that without the 0+ 
in <P{x+Q'^) in the expressions for R{x) and L{x) (i.e. the 
last expressions on the RHS of Eq. (C42)), the pre- 
multiplying the last term on the RHS of Eq. (C75) would 
have been '— ' (compare with Eq. (6.29) in Ref. [28]; see 
also the remarks centred around Eq. (C80) below). This 
follows from the identity (cf. Eq. (C70)) 

2 cos(.t) cos(y) = cos(a; — y) + cos{x + y), 

which would have resulted in +1 instead of the —1 in the 
first term on the RHS of Eq. (C73). 

The Hamiltonian in Eq. (C75) coincides with that in 
Eq. (6.29) of Ref. [28] (taking into account that Jz in the 
latter reference is the equivalent of the V in our present 
considerations [29], that = aJz, and that, for the rea- 
son specified in the previous paragraph, the sign of the 
last term on the RHS of Eq. (C75) is opposite to that of 
the corresponding term in Eq. (6.29) of Ref. [28]). Fur- 
ther, the coefficient of the cosine term in Eq. (29.27) of 
Ref. [4], and that in Eq. (11.14) of Ref. [31], is twice 
too large; this is due to partially neglecting the 2 on the 
LHS of the identity in Eq. (C70) above. In verifying this 



statement, one should note that the ratio V/ {2t) is the 
equivalent of the dimensionless quantity A, Eq. (C104), 
the anisotropy parameter, in Eq. (29.27) of Ref. [4], and 
in Eq. (11.14) of Ref. [31]. 



3. The (exact) bosonized Hamiltonian for small \t'\ 
(for t' = 0) and arbitrary values of t and V 

As we have indicated earlier, the above analysis is valid 
for V sufficiently small with respect to t and t' . Here 
we consider the specific case corresponding to t' = 0, for 
which the model under consideration is identically equiv- 
alent to the XX Z Heiscnberg Hamiltonian for SU(2) spin 
operators. Regarding the specific choice t' = 0, we refer 
the reader to our pertinent remarks in the paragraph 
preceding Eq. (C42) on page 16. Since the transfer ma- 
trix of the zero-field eight-vertex model [42] (defined on 
a two-dimensional lattice) commutes with the Hamilto- 
nian of the one-dimensional XYZ model for SU(2) spin 
operators [43, 44], on the basis of the exact solution of 
the former model [42, 45], one can deduce the low-energy 
properties of the latter model [27, 46]. This enables one 
to deduce the bosonized Hamiltonian of the t-V model 
at half-filling for in principle arbitrary values of V. 
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The bosonizcd Hamiltonian corresponding to the t-V Hamiltonian can be brought into the standard sine-Gordon 
form 



(C76) 



r 



A comparison of this Hamihonian with that in Eq. (C75) 
reveals that for sufficiently small V/t one has 



(C77) 

(C78) 
(C79) 



9 = 


2aV 


(27ra)2 ' 


1 


(-^)"^ 


K ~ 


u ~ 





Later we shall establish (see the discussions following 
Eqs. (C116) and (C123)) that the RHSs of the expres- 
sions in Eqs. (C78) and (C79) are exact to leading order 
in V/t. 

For completeness, we point out that the sign of the 
coupling constant g is not crucial, as it can be reversed 
on effecting the shift transformation: 



21 + 1 



(C80) 



where / is an arbitrary integer (negative, zero and posi- 
tive) . Since the constant term on the RHS of Eq. (C80) is 
a c number, the transformed ^{x) remains to satisfy the 
commutation relations in Eq. (C18); consequently, the 
commutation relation in Eq. (24) remains equally appli- 
cable to the transformed 'P{x) as the untransformed ^(x). 
With reference to Eqs. (C6) and (C42), the transforma- 
tion in Eq. (C80) gives rise to the multiplicative phase 



factors c+*(2'+i)'^/^ and e «(2'+i)'t/4 jj-^ ^i^^ expressions for 
R(x) and L{x) respectively. 

In the following we shall have occasion to refer to the 
Euclidean action corresponding to the Hamiltonian in 
Eq. (C76), which we determine now. With r the imag- 
inary time, Eq. (C45), on the basis of the Hamilton 
equation of the classical mechanics of point particles, i.e. 
(Eq. (40.4) in Ref. [47]) 



dH 



(C81) 



one has 



^ (x, —ir) = idr'P{x, —ir) 



sn 



5n(x, —it) 



(C82) 



where Ti stands for the Hamiltonian as expressed in terms 
of the imaginary-time Heisenberg-picture counterparts of 
TI(x) and (l>(x). From Eq. (C76) one obtains that 



{x, —it) — uK n{x, —ir). 



(CSS) 



For the Euclidean action one thus has (cf. Eqs. (2.1) and 
(40.2) in Ref. [47]) 



S ^ - 



= dx dr 



dr |77(x, -ir) S (x, -ir) - H | 
1 r 1 



{^[^{dr'P{x,-iT)y + ^{^J{x,-^T)y] +9COs{VT6^S{x,-tT))}, (C84) 



where (3 = 1/T (with h = = \). In the second ex- 
pression on the RHS of Eq. (C84), / stands for /^^/(f. 
Introducing 



S{x,-iT) = S{x,~iT), 

n{x,—iT) = n{x,~iT), 



where 



UT, 



(C85) 
(C86) 



(C87) 



and subjecting the fields ^(x, —ir) and fl{x, —ir) to the 
canonical transformations (see Eq. (24)) 

(C88) 
(C89) 

the action in Eq. (C84) can be expressed as 

S ^So + g I d^r cos((3|(a;, -ir)), (C90) 



(P{x,-iT) Vif f (.T, -it), 

n{x,-iT) -^n{x,-iT), 

V K 
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where 



So 



(driix, ~lT)f+{d.A{x, ~lT)f (C91) 



is the Gaussian action, and 



etc. Using the expression for (3 in Eq. (C92), one deduces 
that for 



K > 



(C97) 



9 



(C92) 



The constant (3 is not to be confused with /? = 1/T. 
Above, 



d^r = da; dr = u dx dr u d^r. 



(C93) 



The two-dimensional integrals in Eqs. (C90) and (C91) 
are over the rectangle [0, L] x [0, [3], where 



(3^ up. 



(C94) 



In the following, where no confusion can arise, we shall 
suppress the bars below the symbols t, g and 'P{x, —it). 

With d denoting the scaling dimension of the field 
cos((3^(a;, — it)) (see Eq. (ClOO)), for which one has (see 
Eq. (C102)) 



47r' 



(C95) 



and D the dimension of the space over which the integral 
in Eq. (C90) is carried out (here D = 2), according to a 
well-known theorem one has that for d > D the perturba- 
tion to the Gaussian action Sq in Eq. (C90) is irrelevant, 
for d < Z? it is relevant, and for d = D marginal (see 
pp. 206 and 207 in Ref. [4]). In the case at hand, where 
D — 2, the perturbation to the Gaussian action 5*0 is 
irrelevant when 



(C96) 



the perturbation to the Gaussian action is irrelevant, for 
if < ^ it is relevant, and for K = ^ marginal. Thus, for 
K > a-iid sufficiently small g, perturbation theory ap- 
plies and the essential physics of the sine-Gordon action 
is determined by the Gaussian action Sq. Amongst oth- 
ers, for K > ^ the single-particle excitation spectrum of 
the system is gapless and thus the system is 'critical' and 
the GS 'disordered' (see later, in particular the remarks 
in the paragraph preceding Eq. (C145) on page 26). With 
reference to Eq. (C78), one observes that for K k, ^ one 
must have V/t k, 37r/2 w 4.71, which is too large for 
Eq. (C78) to be accurate. As we shall see later, K ~ ^ 
exactly corresponds to V/t = 2. Thus, with Vc{t,t') de- 
noting the critical value of V for the one-dimensional 
t-t'-V model on a regular lattice, one has (cf. Eq. (C119) 
below) 



Vcit,0)^2t. 



(C98) 



With 



z^T + ix/u and z'^r' + ix'/u, (C99) 



the scaling dimension d of the field cos(|3^(a;, — ir)) is 
defined according to the following expression: 



(cos(|3<^(a;,-ir))cos((3<^(a;',-ir')))o 1 pti for \z-z'\^oo, (ClOO) 

\z — z'p" 

where C is a constant. The average (. . .)o is with respect to the Gaussian action 5*0, whereby the u's in Eq. (C99) 
are to be momentarily identified with v^, Eq. (C46). To calculate d, one employs the identity cos(a;) = |(e" -|-e~'^). 
Using the left-most identification in Eq. (C47), on the basis of the fact that (cf. Eqs. (22.8) and (22.9) in Ref. [4]; see 
also Appendix C in Ref. [28]) 

{e'^'^(^'^^ e»p*(2',2')\^ = for - ^ oo, (ClOl) 

a 

one obtains that 

(cos(|3<^(x,-ir))cos(|3<^(x',-iT')))o = i| (c'P*(---") c-'P*(-''«'))o -f (e-''P*(^'^") e*P*(^''^"))o} 

for - > 1, \z-z'\ > a. (C102) 
a 
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Comparing the last expression in this equation with the RHS of Eq. (ClOO), one arrives at the expression for d in 
Eq. (C95). 



4. Exact results for K, u and M 

Following Eq. (C92), (3 is seen to be determined by 
K, Eq. (C76). For determining K, one employs the cor- 
respondence between the t-V Hamiltonian (or the t-t'- 
V Hamiltonian at half-filling and for small values of t' 
— see our pertinent remarks in the paragraph preced- 
ing Eq. (C42) on page 16) and the XXZ Heisenberg 
Hamiltonian for SU(2) spin operators (also known as 
the Heisenberg- Ising Hamiltonian), which is established 
through the following identifications (compare Eq. (4) 
with Eq. (6.10) in Ref. [28]): 

J.^y 2t, ^ V. (C103) 

Here, Jxy denotes the common value of Jx and Jy. One 
defines 

J, V 

A=-^^- (C104) 

J xy 

which measures the anisotropy of the exchange parameter 
in the z direction in comparison with those in the x and 
y directions. 

For the XXZ Heisenberg Hamiltonian, in the thermo- 
dynamic limit, one can show that for |Z\| < 1 the system 
is 'critical' and the GS 'disordered' (see later, in partic- 
ular the remarks in the paragraph preceding Eq. (C145) 
on page 26); at Z\ = — 1 the GS of the system is an 
isotropic Heisenberg ferromagnet, and at Z\ = -1-1 it is 
an Heisenberg anti- ferromagnet. For A < —\ the GS is 
an Ising ferromagnet (FM), and for Z\ > 1 it is an Ising 
anti- ferromagnet (AFM). Whereas the Ising FMic GS is 
gapless, the Ising AFMic GS is gapped. The Ising AFMic 
GS has a broken translational symmetry, characterised 
by a non-vanishing staggered magnetization the or- 
der parameter of the GS (for the relationship between 9Jl 
and 5n as introduced in Eq. (20), see the remarks follow- 
ing Eq. (C63) above); the continuous symmetry of the 
GSs corresponding to Z\ < 1 (i.e., the chiral symmetry 
of the underlying quantum sine-Gordon Hamiltonian) is 
for A > 1 spontaneously broken and reduced to a dis- 
crete Z2 symmetry. Thus, the unit cell of this broken- 
symmetry GS is twice as large as that of the uniform GSs 
corresponding to Zi < 1 (Ch. 11 in Ref. [31], Ch. 6, in 
particular Fig. 6.1, in Ref. [28] and Ch. 29 in Ref. [4]). 

Below we consider in some detail the cases correspond- 
ing to |Z\| < 1, Z\ < -1 and A>1. 

a. The case of \A\ < 1 
(including a discussion of the XY Z model, with Jx 7^ Jy) 

The value of K is deduced from the leading-order 
asymptotic term of the correlation function 

(5-(x,-*T)^-(0,0)),tag, 



or of that of 

(5+(x,-zr)5-(0,0)), 

in the asymptotic region |z| = |r + ix/u\ —f 00 (cf. 
Eq. (C99)). Here S'{x,—iT) is the z component of 
the SU(2) spin operator in the continuum limit, and 
S^ {x, ^ir) and S~ {x, —it) are respectively the rasing 
and lowering SU(2) spin operators, also in the continuum 
limit. The latter three spin operators can be expressed 
in terms of the field operators <?(x, — ir) and 0{x, —it), 
whereon the above-mentioned asymptotic terms are rel- 
atively straightforwardly calculated (see Ref. [48] and 
Ch. 6 of Ref. [28]; for a complementary calculation re- 
garding the static counterparts of these dynamic corre- 
lation functions, see Refs. [49, 50]); on replacing (. ..) 
by (. . .)o, these asymptotic terms are readily calculated 
along the same lines as those leading to Eq. (C102) above. 

The subscript 'stag', for 'staggered', refers to that part 
of the correlation function (S'^(a;, — ir)S'^(0, 0)) whose 
amplitude is modulated by the phase factor e^'*^^^ (at 
half-filling, for which = Eq. (C3), this phase fac- 
tor gives rise to a sign alteration for the j in a; = ja in- 
creasing in steps of unity) . We note that the asymptotic 
series expansion of {S''{x, — ir)S'^(0, 0)) corresponding to 
[zl — > 00 contains a term that decays like 1/|2|^, so that, 
in view of Eq. (C106) below, the leading-order term in the 
asymptotic series expansion of {S^ {x, —it) S^ {0,0)) stag, 
specific to |z| — > 00, may not be the leading-order term of 
the asymptotic series expansion of {S^{x,—iT)S^{0,0)). 
This is in fact that case for — 1 < Z\ < 0, for which one 
has 1 < K < 00; for < Z\ < 1, in contrast, one has 
i < X < 1 (see Eqs. (C108) and (C116)). 

For the explicit asymptotic expressions of the first (sec- 
ond) of the above-indicated correlation functions, we re- 
fer the reader to Eqs. (16) and (21) [51] (Eqs. (15) and 
(15')) of Ref. [48] (compare with Eqs. (29.37) and (29.38) 
in Ref. [4], and with Eq. (6.38) in Ref. [28]). In this con- 
nection, wc note that 

A^^^, (C105) 

where 6 is the symbol employed in Ref. [48] {9 is denoted 
by 77 in Refs. [49, 50]). For |z| 00 one has 

\{S''{x,-tT)S'{o,o)U.g\ ~ (C106) 



\{S+{x,-tT)S-{0,0))\^j-^^, (C107) 

where Ci and C2 are some finite constants. These power- 
law decays are indicative of the mass-less spectrum of the 
uniform, or 'disordered', GSs corresponding to |Z\| < 1, 
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which, as we briefly indicated in the previous paragraph, 
correspond to K > ^ (cf. Eq. (C97)). This is a pecuhar- 
ity of the XX Z model, since for the XYZ model, the 
deviation of from Jy gives rise to a massive spectrum, 
as we shall discuss later in this section (see Eqs. (C142), 
(C143), (C144) and (C163) below). 
For |Z\| < 1, A is expressed as 



A = +cos(h), < h < tt. 



(C108) 



where \x is one of the elliptic-function parameters that one 
encounters in the exact solution of the 8- vertex model on 
a square lattice in the thermodynamic limit (see Ref. [46] 
and Ch. 10 in Ref. [45]). 

The plus sign on the RHS of Eq. (C108), as opposed 
to the minus sign in the counterpart of this expression in 
Eq. (10.16.8) of Ref. [45], calls for an explanation. The 
8- vertex model on a square lattice is solved in a so-called 
'fundamental region' (FR) of the parameter space (the 
parameters being conventionally denoted by wi, 102, 
and W4 — Eqs. (10.2.16) and (10.1.2) in Ref. [45]) and 
the solution for other regions of this space are deduced 
from the symmetry of the partition function of this model 
under some well-specified transformations of these pa- 
rameters (Eq. (10.2.17) in Ref. [45], and Appendix A in 
Ref. [46] ) . The essence of this symmetry property is en- 
capsulated in a 'rearrangement procedure' which for the 
specific case of the XYZ model is very simple. For this 
model, the FR is defined as the region for which (Sec. 
10.15 in Ref. [45]) 



\Jv\ <Jx< -Jz 



(C109) 



holds. For the cases where the parameters {Jx, JyT-Jz} 
do not satisfy the inequafities in Eq. (C109), one de- 
fines the 'rearranged' parameters {J^,Jy,J|}, satisfy- 
ing these inequalities, obtained through a permutation 
of {Jxt Jyi Jz} in combination with possibly negation of 
a pair of the original parameters. The expressions pre- 
sented in Eq. (10.16.8) of Ref. [45] correspond to 

{Jx,Jy,Jz)^{l,hA), (CllO) 

for which in the case of |Z\| < 1 one has: 



(Ji,j^,JD = (i,-A-i). 



(Clll) 



This in turn leads to il ~ —A and A^ = —1 (for the 
definitions of il and A^, sec Eq. (10.15.5) in Ref. [45]). 
In our case, in contrast, one has (upon an appropriate 
normalisation of the parameters - see Eqs. (C103) and 
(C104); compare the Hamiltonian in Eq. (19.1) of Ref. [4], 
or that in Eq. (6.2) of Ref. [28], with the Hamiltonian 
in Eq. (10.14.1) of Ref. [45] and in doing so note that 
5^ = ifff , where v = x,y,z): 

{Jx,Jy,Jz) = (-1,-1, -zl), 

which in the case of |Z\| < 1 leads to 



(C112) 



In our case we thus have /I = +A and A^ = —I (compare 
the latter with the one presented following Eq. (Clll) 
above). 

For completeness, in the case of |Z\| < 1, for the el- 
liptic modulus K (often denoted by k) one has k = 1 
(Eq. (10.15.6a) in Ref. [45]). Due to the fact that 

snh(A, k) = —i sn(iA, k), 

which for K = 1 is identical to tan(A), from Eq. (10.15.6b) 
in Ref. [45] one immediately arrives at Eq. (C108). In 
doing so one makes use of Eq. (10.15.12) in Ref. [45], in 
which /' = K{k'), where k' = (1 — k^)^/^ is the elliptic 
conjugate modulus; for k = 1, k' = 0, so that /' = tt/2 
(whereby A = \x/2). 

Above, sn(u, k) is the Jacobian elliptic sinam function 
(Sec. 22.11 in Ref. [37]) and K is the complete elliptic 
integral of the first kind (Ch. 15 in Ref. [45]), or "the 
constant K" [K = K{k,)) corresponding to k (Sec. 22.3 
in Ref. [37]); by definition, sn(if , k) 1 (Sec. 22.3 in 
Ref. [37]). To avoid confusion with the K as encountered 
in e.g. Eq. (C76), in what follows we shall adopt the 
notation in Ref. [4-5] and denote K = K{k) by I and 
K' = K{k') byl'. 

As we shall have occasion to refer to some results in 
Ref. [23] regarding elliptic functions, we point out that 
whereas in our notation is the parameter of ellip- 
tic functions and k'^ the complementary parameter, in 
Ref. [23] (in particular in Chaps. 16 and 17 herein), m 
is the parameter of these functions and mi = 1 — m the 
complementary parameter (Sec. 16.1 in Ref. [23]). 

For later use, we further remark that in the theory 
of elliptic integrals, one encounters the nome q = q{m) 
and complementary nome qi = q{mi), defined accord- 
ing to (items 17.3.17 and 17.3.18 in Ref. [23]; see also 
Eq. (15.7.6) in Ref. [45]) 

q = q{m) = exp[— 7r/'//], qi = q{mi) = cxp[— tt///']. 

(C114) 

For q one has the following expansion (item 17.3.21 in 
Ref. [23]): 



'^("^)=16+Hl6 



for jmj < 1. 



-«HT6. 

(C115) 

This series is very useful in deducing the asymptotic se- 
ries expansion of q for m —^ and that of qi for to — )■ 1 . 

The elliptic-function parameter |i. in Eq. (C108) is re- 
lated to K (as encountered in Eq. (C76)) as follows [48] 
(cf. Eqs. (C105) and (C162)): 



7r(l 



1 

2K 



K 



2(71 - h) 



(C116) 



{Jl,JlJl) = {i,A-i). 



Thus the K corresponding to the cases where |Z\| < 1 
is deduced by obtaining the value of |j, from Eq. (C108) 
and substituting the resulting value into the right-most 
expression in Eq. (C116). From Eq. (C116) one ob- 
serves that the non-interacting case, for which K ^ \, 



(C113) Eq. (C78), corresponds to \i 



and. in view of 
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Eq. (C108), to Z\ = 0. Following Eq. (C104), indeed 
Zi = corresponds to V — 0. 

Making use of the series expansion 



/ \ 1 3 

arccos X = x x 

^ ' 2 6 



. for X — > 0, 
from Eqs. (C108) and (C116) one deduces that 



2A 4Z\2 
A = 1 \ — 



(C117) 



to be compared with the result obtained from the expres- 
sion in Eq. (C78), i.e. 



K = 1- 



2A 6Z\2 



(C118) 



One observes that to the linear order in A ^V/ (2t) the 
two expressions coincide. 

From Eqs. (C116), (C108) and (C104) one deduces that 



K 







1 



V = 2t. (C119) 



In particular, the latter exact equality (cf. Eq. (C98)) 
deviates from the perturbative result V ~ (obtained 
from Eq. (C78)) to which we referred in the text following 
Eq. (C97) above. 

To make contact with Ref. [48], we note that the 
expression in Eq. (C116) is the equivalent of that in 
Eq. (21') of Ref. [48] where J^y is identified with unity 
and thus where Jz is the equivalent of the A in our 
present considerations (cf. Eq. (C104)); for clarity, one 
should recall Eq. (C105) above and note that arcsin(Z\) + 
arccos(Z\) = 7r/2. To make contact with Ref. [28], we 
note that the counterpart of the parameter |x in our 
present considerations is the 7r(l — 0^) (and not tt/?^) 
in Ref. [28] (see Eq. (6.31) herein), where the latter (3 is 
not to be confused with the [3 in our present considera- 
tions, which stands for 1/T; the 0^ of Ref. [28] is thus 
equivalent to the 9 in Ref. [48] ; this is easiest appreciated 
by comparing the expression 1/K = 2/3^ in Eq. (6.31) of 
Ref. [28] with Eq. (C105) above. 

Having determined the functional form of the exact K 
corresponding to |Z\| < 1, we now proceed with doing the 
same for the exact u, for which for small values of V/t 
one has the expression in Eq. (C79). 

For the excitation energy AE corresponding to the case 
|Z\| < 1 one has [46] (see Eq. (7.10b) herein) 



^ ^ TT sin(Li) /, . , , , 
AE = t (I sin(agi)| 



sm 



(ag2)|), (C120) 



where \i is the parameter deduced from Eq. (C108) and 
a the lattice constant (identified with unity in Ref. [46]). 
We note that the last expression in Eq. (7.10b) of Ref. [46] 
contains typing errors: instead of sin (71 and sin 52, one 
must have | sin qi \ and | sin \ . This fact can be verified 
by realising that for |Z\| < 1, the ki and fc2 in Eq. (7.8) 
of Ref. [46] are equal to 1 whereby the terms involving 



cos^ qi and cos^ q2 reduce to | sin qi \ and | sin 52 1 respec- 
tively. Physically, this fact can be appreciated by the 
realisation that the excitation energy AE cannot be neg- 
ative (see the discussions concerning elementary excita- 
tions in metals in Sec. 1.1, p. 10, of Ref. [52]). 

For clarity, comparing the Hamiltonian in Eq. (19.1) 
of Ref. [4] (or that in Eq. (6.2) of Ref. [28]), with that 
in Eq. (7.1) of Ref. [46], and in view of S'f = iaf , v = 
x, y, z, we point out that the — in Eq. (7.8) of Ref. [46] 
(— is identified with 1 in Eq. (7.10b) of this reference) is 
to be replaced by t and not 2t (cf. Eq. (C103) above); on 
the absolute scale, the magnitudes of energies in Ref. [46] 
are twice as large as those in our present considerations. 
This clarifies the reason for the t, as opposed to 2t, on the 
RHS of Eq. (C120) (compare with Eq. (6.31) of Ref. [28] 
where one also encounters ^Jxy, instead of Jxy, in the 
expression for u). 

With 



q = qi+q2 



(C121) 



denoting the wave number of the low-lying elementary 
excitations in the system, from Eq. (C120) one obtains 
that 

d . ^ TT Sin(|j.) / . / r 1N\ / r i\ 

■7T-AE = at sgn(sm(a g - 92 )) cos(a g - 92 )• 

oq |X ^ ' 

(C122) 

One can readily verify that the minimum of AE corre- 
sponds to qi ~ q2 = rmr, where m ~ 0,±1,... . For 
q = 2m7r -|- 0+ the expression in Eq. (C122) thus yields 
the value for u. One thus obtains that (cf. Eq. (29.31) 
in Ref. [4] and Eq. (6.31) in Ref. [28]) 



TT sin(|j.) 



(C123) 



where we have used Eq. (C40). For small values of Z\ = 
V/{2t), from Eqs. (C108) and (C123) one obtains that 



u ^ Vf(1 



2A (^2 _ 8)^2 
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which is to be compared with 



Wf(1 + 



2A 2Z\2 



(C124) 



(C125) 



which is obtained from the expansion of the expression in 
Eq. (C79) in powers of A. Thus indeed to linear order in 
V/t the expression in Eq. (C79) coincides with the exact 
expression in Eq. (C123). 

From Eqs. (C119) and (C123) one deduces that 



K 



(C126) 



This result deviates from the perturbative result u w 2vf 
that one obtains from Eq. (C79) and V w ^t, the latter 
being deduced from Eq. (C78) and the condition K = ^. 
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The value of the gap in the singlc-particlc excitation 
spectrum corresponding to Z\ > 1 can be directly de- 
duced either from the expression for the inverse corre- 
lation length or from that for the singular part of 
the Gibbs, or Helmholtz, free energy near criticality (see 
Eqs. (10.12.22a) - (10.12.23b) and (10.14.39) - (10.14.43) 
in [45]). In general, the energy gap is proportional to 
where z is the dynamic scaling exponent (Eq. (1.3) 
in Ref. [20]). In our case, the Lorentz invariance of the 
theory under consideration (see Eq. (C39)) implies that 
z is equal to unity (p. 88 in Ref. [20]) so that the en- 
ergy gap is directly proportional to For a general 
discussion of the inverse correlation length, the reader is 
referred to Sec. 16.2 in Ref. [53] and Sec. 5.2 in Ref. [54]. 

In anticipation of the calculation of the mass gap Ai in 
the broken-symmetry GS corresponding to Z\ > 1 in Sec- 
tions C 4 c and C 5, below we establish the critical nature 
of the GS of the system imder consideration for \A\ < 1 
(a fact that is apparent from the power-law decays of the 
correlation functions in Eqs. (C106) and (C107) above) 
by considering the case Jx = Jy as the limit of Jx i Jy 
To this end, we consider the case corresponding to 



\Jz\ <Jy< Jx- 



(C127) 



Following our earlier considerations (regarding the plus 
sign on the RHS of Eq. (C108)), one has 

(J^,J^,JD = (J„-J.,-J.), (C128) 
so that (Eq. (10.15.5) in Ref. [45]) 

r, = -^, A = (C129) 

Jy Jy 

In view of Eq. (C127), it will be convenient to introduce 

e = Jx-Jy (C130) 

and consider the case e j 0. For convenience, in what 
follows we assume that 



''/Jy 



(C131) 



Following Eq. (10.15.6a) in Ref. [45] and making use 
of the expressions in Eq. (C129) for Fi- and A^, one has 



2kV2 

1 + K 



= 1 -7^ + 0(67?) +0(77^), 



(C132) 



from which one obtains that (k is bound to satisfy < 
K < 1, Eq. (10.15.7) in Ref. [45]) 

K = l-2(2r/)i/2 + 0(?7). (C133) 

Consequently, 

k' = {1- K^ fl^ = 2 (2r,)i/4 + 0(7,3/4). (C134) 



Making use of 

snh(A,K) = tan(A) + |-Asec2(A)+tan(A)| {2r]f^'^^0{T]), 

(C135) 

and expressing A as 



A = Ao + (5A, 



(C136) 



where Aq is the solution of Eq. (10.15.6b) in Ref. [45] 
corresponding to e = (and thus to 77 = 0) — that is, Aq 
is the solution of 



tan(Ao) 



1 + 



■ 1 J z I Jy 

after some algebra for 5\ one obtains that 



cos(2Ao) 

(C137) 



<5A = Ao (277)1/2 + 0(7,). 



(C138) 



We note that < A < /' (Eq. (10.15.7) in Ref. [45]), 
where, in view of Eq. (C134), for /' one has the expression 
in Eq. (C147) below. Thus, < Ao < 7r/2, which is in 
conformity with the first equation in Eq. (C137) of which 
the RHS is positive. 

In order to establish the critical nature of the sys- 
tem under consideration for Jx = Jy and |Z\| < 1, we 
need to calculate the quantity t, defined in Eq. (10.12.2) 
of Ref. [45], which is vanishing on the critical surface; 
in general, this quantity vanishes linearly with T — Tc, 
where is the critical temperature; (T — Ts)/t is posi- 
tive (Sec. 10.12 in Ref. [45]). In this connection, we note 
that the parameters wi, W2, and W4, or a, 6, c and 
d, of the 8-vertex model are functions of temperature 
T, Eqs. (10.1.2) and (10.2.1) in Ref. [45]. According to 
Eq. (10.12.3) of Ref. [45], one has 



K + = 2-4t. 



(C139) 



Calculation of t from this expression is mediated by 
the calculation of the quantity p for which one has 
(Eq. (10.12.5) in Ref. [45]) 



(1-At2)2 (1-^2)3 



0((1 



162 162 

(C140) 

where we have employed the series expansion in 
Eq. (C115) (for qi see Eq. (C114)). Expressing k as 
K = kq+Sk, with kq denoting the solution of the equation 
(1 - ^2)2 ^ lQ2p^ 

no = il- 16pi/2)i/2 ^ 1 _ 8pi/2 _ + 0(^^/2), 

from Eq. (C140) one obtains that 

K= l-8pi/2-|-32p+0(p3/2). (C141) 

Substituting this expression into Eq. (C139), for k 
1, and thus p — > 0, Eq. (C140), one deduces that (cf. 
Eq. (10.12.19) in Ref. [45]) 



t 



-16p+0{p^) 



P 



1 

16 



t + 0(t2). (C142) 
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Making use of Eqs. (C140), (C133) and (C130), one has 
(cf. Eq. (10.14.37) in Rcf. [45]; see also Sec. 4.5 in 
Rcf. [32]) 



^ '^^ + 0{if/')+0{er,/Jy). (C143) 



16 J2 - J2 



For the inverse correlation length ^ ^ near criticality, 
one has (Eq. (10.12.23a) in Ref. [45]) 



'1 ^ £(_t)^/2^^ 



(C144) 



where C > is a finite dimensionless constant and a the 
lattice constant. The following four remarks arc in order: 

First, since Eq. (C144) applies for t < 0, and since 
(T — Tc)/t is positive (see the previous paragraph), it 
follows that the case at hand, which is characterised by 
Jx > Jy > \Jz\, corresponds to the 'ordered region' (T < 
Tc) of the 8- vertex model. 

Second, the mathematical details underlying the ex- 
pression in Eq. (C144) are strictly valid for A < ^I' (see 
pp. 241 and 242 in Ref. [45]). Baxter however has shown 
that Eq. (C144) itself is valid for A < §/', i.e. for h < ^ 
(see pp. 240 and 253 in Ref. [45]). Nonetheless, Johnson 
et al. [46] have demonstrated that Eq. (C144) gives the 
correct critical behaviour even for |j, > ^ (see p. 253 in 
Rcf. [45]). 

Third, Eq. (C144) is not vahd for |j. = 0; in such 
case, vanishes exponentially as criticality is ap- 

proached; the appropriate expression for in such 
case is that presented in Eq. (8.11.24) of Ref. [45] (see 
p. 272 in Ref. [45]). We shall encounter this behaviour in 
Eq. (C187) below. The renormalization-group analysis 
in Sec. C 5 reveals that this behaviour is specific to the 
states inside the cross-over regime (CO), adjacent to the 
strong-couphng (SC) regime, to be discussed in Sec. C 5 c 
below. In contrast, the behaviour of as described by 
Eq. (C144) is specific to the states deep inside the SC 
regime, to be discussed in Sec. C 5 a. 

Fourth, Eq. (C144) apphes strictly for ■k/{2\i) different 
from an integer; for tt/{2[x.) an integer, the expression 
for acquires a logarithmic correction (see p. 252 in 
Ref. [45]; see also Sec. 8 in Ref. [42]). 

One observes that for i Jy, one has p i (see the 
inequalities in Eq. (C127)), whereby t t 0, Eq. (C142). In 
the limit = Jy, while | < J^ = Jy (cf. Eq. (C127)), 
or |Z\| < 1 in the hght of Eq. (C103), t is vanishing, sig- 
nifying the system under consideration as being critical. 
For Jx > Jy > \ Jz\, P non- vanishing and unless |j, = 
(see the summary in the paragraph following Eq. (C158) 
on page 27), is non- vanishing so that the spectrum of 
the system is gapped (see the remarks in the paragraph 
following Eq. (C126) on page 24). 

We proceed now with the calculation of \i, for which 
one has (Eq. (10.12.5) in Ref. [45]) 



(C145) 



For Jx = Jy and | Z\ | < 1 we should recover the expression 
in Eq. (C116) above. Making use of (item 17.3.11 of 
Rcf. [23]) 

I{rn) = ||l+ (i)'m+ (i-|)'m2 + ...}, |m| < 1, 

(C146) 

from Eq. (C134) (and mi = k'^) one obtains 

/' = | + |(2r;)V2 + 0(^). (C147) 

Thus, with 

Ho - ^ Ao = 2Ao - 2Ao(277)i/2 + 0(,y) , (C148) 
and (cf. Eq. (C138)) 

<5h= ^<5A = 2Ao (277)1/2 + o(^), (C149) 



one has 



(C150) 



so that on account of Eq. (C137) (recall that < Aq < 
V2) 

(C151) 



2y. 4Ao + 0(?7) 2 cos-i(- Jy) + 0(77) ' 

For Jx = Jy (= Jxy, Eq. (C104)), i.e. for e 0, which, 
following Eq. (C131), imphes = 0, Eq. (C151) yields the 
exact result A = -cos(^), Eq. (10.16.8) in Ref. [45]; the 
origin of the minus sign in this expression, which devi- 
ates from the plus sign in Eq. (C108), lies in the fact that 
in arriving at Eq. (C151) we have been dealing directly 
with [Jx, Jy,Jz), ordered according to Eq. (C127), lead- 
ing to Eq. (C128), where one has the signature {-+, — , — ). 
This is to be contrasted with the signature of the corre- 
sponding parameters in Eq. (C113), which is (+,+,—). 
For further relevant details, the reader is referred to the 
discussions in the paragraph following Eq. (C108). 
For the cases where 



Jz 



one has 



1 - 7 with < 7 < 1, 



A.^(?)"Vo(7"') 



(C152) 



(C153) 



It would therefore appear that for 7 | 0, 7r/(2|j.) would 
approach a large value, of the order of l/ry, which is in- 
dependent of 7. This is however not the case, as 77 is a 
non-trivial function of 7; from Eqs. (C131) and (C152) 
one obtains that 



2Jy 7 



(C154) 



which implies that as 7 approaches e from above, the 
condition rj I, Eq. (C131), which underlies our above 
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considerations, will eventually fail to be valid. It follows 
that the validity of Eq. (C151) is conditioned on e ^ 7. 
Thus, although 7 can approach zero, it my not come 
closer to zero than e = Jx — Jy in order for Eq. (C151) 
to be applicable. 
For the cases where 



one has 



Jz 

J il 



Ao = ^ - 4 



-1+7 < 7 < 1, 



'y\ 1/2 

+0(7^/2), 



so that 



^ = Ho + '5H = 7r- (27)1/2 



from which one obtains that 
TT 1 1 



2h 



Ztt 



(C155) 



(C156) 



(C157) 



(C158) 



It can be shown that Eq. (C154) applies also here, so that 
the smallness of rj similarly requires e <C 7. 

Summarising, since the exponent 7r/(2|j,) as deduced 
from Eq. (C151) is finite, with reference to our remarks 
in the paragraph preceding Eq. (C145) on page 26, we 
conclude that indeed for > Jy the inverse correlation 
length is non- vanishing and thus the excitation spectrum 
of the system is gapped. In contrast, for = Jy (i.e. 
for e = 0) the system is critical (t = 0) and the GS 
is 'disordered' (^"^ vanishing). This fact is reflected in 
the power-law decays of the GS correlation functions in 
Eqs. (C106) and (C107). 

For Jx > Jy (generally, for Jx ^ Jy) and j^j — > 00 one 
has [48] (see p. 3912 herein) 



{S+{x,-iT)S- {0,0)) 



C3 _ C3 

= J^l/(2K)- 



(C159) 



instead of the expression in Eq. (C107). Here C3 is a 
finite constant and M the correlated mass gap, for which 
one has [48] 



X cx TW^, 
in which (cf. Eq. (C130)) 

= i(J. - Jy) = |, 



(C160) 



(C161) 



the basal-plane anisotropy, is the uncorrelated mass gap. 
One has the following scaling law [48] (see Table I herein) 



K 



2(1 - 61) 2m. 2K - 1 



(C162) 



To appreciate the expressions in Eqs. (C160) and (C161), 
as well as the equivalence of v with 7r/(2|j.) (the second 



equality in Eq. (C162)), one should realise that on ac- 
count of (as regards z = 1, see the remark in the para- 
graph following Eq. (C126) on page 24) 



(C163) 



the expression in Eq. (C160) represents the same fact as 
expressed by Eq. (C144). In this connection, note that 
-t = 2?7 + 0{rf/^), Eqs. (C142), (C143), and that 7? (x e, 
Eq. (C131), so that — t oc A^o near criticality. With 
reference to Eq. (C158), we note that since for e <C 7 
and 7 i one has v [ ^, li follows that for e ^ 7 and 
7 j the constant K approaches -l-oo from below; this is 
evident from v ^ 5 + 31^7 Eq. (C162), which applies for 
\ 00 and which implies i' > ^ for K positive and 
large. 

In the case of the XY model, for which = 0, one 
has the mean-field result ly ^ 1 [48, 55, 56] (see also 
page 270 in Ref. [45]), which, following to Eq. (C162), 
correctly corresponds to K = 1, \x = ^ and 6 — ^. The 
result 1^ = 1 immediately follows from Eq. (C151) for 
the isotropic XY model [55], for which e = 0, and thus 
7] = 0. For the anisotropic XY model [56], Eq. (C151) 
would incorrectly suggest that v deviated from 1. One 
can explicitly demonstrate that for Jz = 0, whereby F,- = 
(see Eq. (C129)), Eq. (10.15.6b) in Ref. [45] yields the 
exact result: 



(C164) 



from which and Eq. (C145) one immediately obtains = 
1. The result in Eq. (C164) follows from snh(A, k) = 
— i sn{iX, k) and the exact result {Example 2 in Sec. 22.41 
of Ref. [37]) 



sn(ii/', k) = IK ^1"^ . 



(C165) 



We note in passing that, on account of r'r = for 
Jz = 0, one has further the following exact solution of 
Eq. (10.15.6a) in Ref. [45] (cf. Eq. (C133)): 

K = 2Z\2-2|A-|(/i?-l)'/'-l, 

where A-c is specified in Eq. (C129). Following the in- 
equalities in Eq. (C127), one has Z\r < ^ 1 so that indeed 
< At < 1, Eq. (10.15.7) in Ref. [45]. 

Below we consider the cases corresponding to |Z\| > 1. 
In doing so, we deal separately with the cases correspond- 
ing to Z\ < —1 and A > \. In what follows, we explic- 
itly consider the case corresponding to the expression in 
Eq. (C112) above. In view of the definition of the FR in 
Eq. (C109), in the following we shall view the condition 
Jx = \Jy\^^ representing 



V 



u:. 



0^ 



(C166) 
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6. The case of A < ~1 

Following Eq. (C112), in this case one has 

(J^,J,^J^) = (1,-1,Z\), so that r,^~l, A = A. 

(C167) 

Since r, = -1, from Eq. (10.15.6a) in Rcf. [45] one ob- 
serves that in the present case 



K = 0. 



(C168) 



With reference to Eq. (10.15.6b) in Ref. [45], which in- 
volves we should consider k = as the limit of 
K I 0. In this hmit, from Eq. (10.15.6a) in Rcf. [45] one 
has 



1 i-r2 



4 z\2 - n " 



(C169) 



Making use of snh(A, k) ~ sinh(A) for k | 0, from 
Eq. (C169) and Eq. (10.15.6b) in Rcf. [45] one obtains 
that 



sinh(A) 



2(Z\2-r2)l/2 



(C170) 



1 + r, 

According to this expression, and in view of Eq. (C166), 

A = +oo. (C171) 

Since n i 0, for the nome q{m) [m = k^) to leading 
order one has q{m) ^ Eq. (C115), on account of 

Eq. (C114) one deduces that 



I TT \kJ 



(C172) 



From this, / ~ tt/2 for k j (item 17.3.22 in Rcf. [23]), 
and Eq. (C146) one thus has 



ttA 



ln(4/K) 



On account of sinh(A) ~ | for A 
Eq. (C170) one deduces that 

'4(Zi?- 1^2)1/2 



A - In 



i + r. 



(C173) 
oo, from 

(C174) 



Conform Eq. (C171), one observes that A indeed diverges 
(logarithmically) for I —1. From the expression in 
Eq. (C174) and that in Eq. (C169) we thus arrive at 



A - In ( - 



according to which and Eq. (C173) (sec Eq. (10.16.8) of 
Ref. [45]) 



lim |x = TT. 



(C176) 



Since for \x = tt the expression in Eq. (C108), which is 
specific to all A € [—1,1], yields A = —1, it follows 
that |J, is a continuous function of A across the boundary 
A = —1 whose derivative with respect to A is however 
finitely discontinuous at this boundary. 

The phase transition at Z\ = — 1 is a first-order one, 
discussed in Sec. 8.11 of Ref. [45]. In Ch. 8 of this refer- 
ence, ice-type models are considered and the phase rel- 
evant to Z\ < — 1 (that is, zi > 1 in the convention of 
Ref. [45]; see our remarks concerning the plus sign on 
the RHS of Eq. (C108)) is the ferroelectric one (i.e. the 
Ising FMic state in the magnetic language — see the sec- 
ond paragraph of Sec. C4, on page 22), corresponding 
to ei < £3, £5 (pp. 156 and 157 in Ref. [45]); reahs- 
ing that the Boltzmann weight d of the 8-vertex model 
is vanishing in the case corresponding to Jx = Jy (the 
Hcisenberg- Ising chain; sec pp. 271 and 272 in Rcf. [45]), 
one observes that indeed for d = the ferroelectric phase 
of the 8-vertex model (phase / as specified on p. 246 
of Ref. [45])) coincides with that of ice- type models; in 
this phase, for the Boltzmann weights a, b and c one has 
a > b + c so that the phase boundary of the ferroelectric 
phase is determined by the condition a = b + c. In this 
case, T - Tc is proportional to (Eq. (8.11.4) in Rcf. [45]) 



t = 



b + c — a 



(C177) 



and for the free energy per site / one has (p. 157 in 
Rcf. [45]) 

r ei-ifcBrct + 0(fV2), t>0, 
/ = < (C178) 
[ei, t<0, 

establishing that indeed df /dt is discontinuous at t = 0, 
or A = —1. Approaching the phase boundary from the 
region t > 0, from Eq. (C178) one further observes that 
the specific heat diverges like t~^/^. 

Since in approaching the ferroelectric phase (our Ising 
FMic phase) one has (Sec. 8.11 in Ref. [45]) 



^ I TT and w I for T 4 Tc 



(C179) 



where w is an elliptic-function parameter, on writing 
(Eq. (8.11.2) in Ref. [45]) 

[i=n-5, w = -n + e, (C180) 

for the Boltzmann weights a, b and c one has 

a : b : c ^ sm{^[e + 6]) : sm{^[e - S]) : sin((5), (C181) 



(C175) resulting in (Eq. (8.11.6) in Ref. [45]) 



t^_(e_5)J for TIT,. 



(C182) 



It follows that for t | 0, both S and e vanish like t^/2^ jt 
should be noted that the approach p. t tt in Eq. (C179) 
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coincides with that required of |x by the expression in 
Eq. (C108) in order for A to decrease towards —1. 

In the ordered phase, corresponding to t < 0, or 
Z\ < — 1, the correlation length ^ is vanishing, and in the 
ferroelectric phase, corresponding to t > 0, or |Z\| < 1, it 
is infinite (p. 157 in Ref. [45]). The latter observation is in 
conformity with the one arrived at in the previous section 
(see the summary in the paragraph following Eq. (C158) 
on page 27). 



that for the case at hand, where A > 1, A is positive. 
Using the expansion 



cosh(A) - 1 + ^A^ for A 0, 
from Eq. (C185) one thus obtains that 



(C188) 



A - +v/2(zi- 1) for Ail, (C189) 
whereby, following Eqs. (C163) and (C187), 



c. The case of A > 1 

Following Eq. (C112), in this case one has 

(j;, JJ) = (1, 1, ~A), so that r, = 1, A, = -A. 

(C183) 

As in the previous case, here also one has k = 0. Making 
use of Eq. (C170), for the present case we thus deduce 
the following equivalent of Eq. (10.15.6b) in Ref. [45]: 



sinh(A) = (Zi^ - 1)1/2, 



(C184) 



From this and the fact that in the case under considera- 
tion A is positive, one arrives at 



A = +cosh(A). 



(C185) 



This result deviates from its counterpart in Eq. (10.16.8) 
of Ref. [45] by a minus sign on the RHS, which devia- 
tion is due to the same mechanism as we discussed in 
connection with the plus sign on the RHS of Eq. (C108) 
above. 

The asymptotic expression in Eq. (C173) applies 
equally for the present case. Since, however, Eq. (C185) 
implies A to be finite for finite values of A (cf. 
Eq. (C171)), in the present case for any finite value of 
A one has 



lim Lx 



0. 



(C186) 



Since for = the expression in Eq. (C108), which is 
specific to all A G [—1,1], yields A = +1, it follows 
that 10. is a continuous function of A across the boundary 
A = +1 whose derivative with respect to A is however 
finitely discontinuous at this boundary. 

As wc have indicated following Eq. (C144) above, the 
expression for in this equation is invalid for |J. = 0. 
For = 0, and T near Tc, one instead has (Eq. (8.11.24) 
in Ref. [45]; see also Eq. (6.12) in Ref. [46]) 



4 TT^ 



(C187) 



As before, the mass gap is proportional to this expres- 
sion, Eq. (C163). The expression in Eq. (C187) can be 
simphfied in the case of Z\ — 1 <C 1. From Eq. (C185) one 
observes that A approaches zero for All. Further, since 
the RHS of Eq. (C184) is positive for Z\ > 1, one infers 



M (X £ 1 ~ - exp , , 

a 2^2{A~1) 



for A I 1. 



(C190) 

This expression is to be compared with that presented in 
Eq. (C226) below. 



5. The mass gap A4: a perturbative 
Renormalization-Group approach 

Comparing the action S in Eq. (C90) above with that 
in Eq. (10.1) in Ref. [31]), for the following considerations 
we introduce the dimensionlcss coupling constant 



g a 



aV 
27r2u' 



(C191) 



where g, the scaled coupling constant, is defined in 
Eq. (C92) above. Assuming that z is small, one can anal- 
yse the consequences of the cosine term in S by means 
of the perturbative renormalization-group (RG) analysis. 
This analysis is standard and its details can be found in 
e.g. Refs. [4, 28, 31]. The considerations in this section 
follow closely those in Ch. 10 of Ref. [31], and Sec. 2.3.2 
of Ref. [28]. 

As we indicated in Sec. C3, the value of the scaling 
dimension d, Eq. (C95), of the cosine field in Eq. (C90) 
determines the significance or otherwise of the contribu- 
tion to S of the term pre- multiplied by g in Eq. (C90): 
for d < 2 the latter term is relevant, for d > 2 it is irrel- 
evant, and for d 2 it is marginal. The boundary value 
d ~ 2 is more generally expressed as d = D, where D is 
the dimension of the space-time domain over which the 
S under consideration is defined. 

With yl ~ 1/a the actual cut-off of the problem in the 
momentum space (see Eq. (C38) et seq.), one writes 



<?(a;) = <Pa{x) 



\k\<A' 

= ^A'{x)+h{x) 



(C192) 



where yl' < yl so that h{x) amounts to the spatially 'fast' 
mode of '1>a{x). Integrating out the contribution of hix) 
to the action S, Eq. (C90), and expanding the resulting 
effective action to second order in g, with 



dl 



A- A' 
A 



(C193) 
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for 



(C194) 



one obtains the following coupled Kosterlitz-Thouless- 
type RG differential equations: 



dz||(0 

dl 
dz^(/) 

dl 



-4(0, 

-Z||(Oz_l(0- 



(C195) 
(C196) 



In Eq. (C194), A is a well-defined constant whose value 
is of no immediate relevance to our present considera- 
tions (see however the last paragraph of this section, on 
page 32), and d{l) is a function of I which is constrained 
by the 'initial' condition 



d(0) 



(C197) 



where d is the scaling dimension of the cosine field under 
consideration, Eq. (C95). For the 'initial' conditions of 
Eqs. (C195) and (C196) one has 

Z||(0) = zjj = d- 2, z^{0) = z1 ^ VsAz. (C198) 

Note that z(0) = z, Eq. (C191). 

Multiplying both sides of Eq. (C195) by Z||(Z), and 
those of Eq. (C196) by z±{l), subtracting the resulting 
equations, one deduces that 



(0-zi(0 



(C199) 



is independent of /; the quantity is therefore com- 
monly referred to as the 'constant of motion' or the 'scal- 
ing invariant' (this fa is not to be confused with ^ and |j, 
that we encounter elsewhere in this paper). One there- 
fore has 



^ 



\2 



(C200) 



The constancy of fi^, for given values of zjj and z^, 
implies that under the RG transformation the locus of 
Z||(^)e|| + z±{l) ej_ on the Z||-zx plane, with {e||,ex} a 
Cartesian basis of this plane, is a parabola. Owing to 
this fact, under the RG transformation, the trajectory 
Z||(/)e|| + zj_{l) e± associated with (z|j,+z']_) is the ex- 
act mirror image of that associated with (zjj,— z*}^). For 
convenience one therefore only considers the trajectories 
in the upper-half of the Z||-zj_ plane. Consequently, un- 
less we indicate otherwise, in the following zj_(l) > 0, 
(with reference to Eqs. (C198) and (C191), in the fol- 
lowing we are therefore tacitly assuming that ^ > 0). 
In this upper half-plane, the region z±_ < Z|| corresponds 
to the weak-coupling (WC) regime, the region z± > |z||| 
to cross-over (CO) regime, and z± < — Z|| to the strong- 
coupHng (SC) regime [28, 31]. The WC and SC sectors 
are characterised by > 0, that is by real fx, and the 
CO sector by — (a-^ = > 0, that is by purely imaginary 



One readily verifies that the following expressions ex- 
actly satisfy the equations in Eqs. (C195) and (C196): 



Z||(0 
z^{l) 



(icoth((i/ 
(icsch(faZ - 



(C201) 
(C202) 



where is a constant, independent of I. One trivially 
verifies that these expressions indeed yield z|(/) — z^(Z) = 

^2 for aU I (cf. Eq. (C199)). Dividing the first of the 
above expressions by the second, identifying I with 0, 
and making use of the initial conditions in Eq. (C198), 
one obtains that 



cosh 



cosh- f 4=1). 



(C203) 



The gap £g in the single-particle excitation spectrum 
of H close to the CDW transition is obtained as follows. 
Let the gap energy corresponding to / 7^ be denoted 
by £g{1) so that for the actual gap one has £g ~ ^g(O). 
Since £qQ) has the dimension of energy (or inverse-time, 
whereby, and by the Lorentz invariance of our continuum 
theory, it has the scaling dimension z = 1), under the RG 
transformation it transforms according to 



£c.{l)=£cM e' 



£c.=£c.{l) e 



-I 



(C204) 



Following the common practice, we consider the follow- 
ing three cases corresponding to gapped states (below, as 
elsewhere, z± > 0): 

(1) z± <C — Z||, specific to states deep inside the SC 
regime, 

(2) z± = — Z|| , corresponding to states on the boundary 
separating the CO and SC regimes, and 

(3) ^± ^ — 2|h signifying states inside the CO regime 
adjacent to the SC regime. 

Before dealing with the above three cases, we de- 
termine the value for the £g{1) specific to the strong- 
coupling regime, where g, the coupling-constant of the 

non-Gaussian contribution to the S in Eq. (C90), is large 
and P, Eq. (C92), is finite. In this connection, we note 
that g (X 1/it, Eq. (C92), so that large values of g can 
be achieved not only by large values of V, but also by 
small values of u. In view of Eq. (C123), u approaches 
zero according to u ~ ^VpiT^ ~ l^) for M- T How- 
ever, following Eqs. (C92) and (C116), |3 - 27r(up/u)i/2 
for |a, I TT, so that for u approaching zero, the mag- 
nitude of gu^^'^ = 5m(3^ (see Eq. (C210) below) only 
weakly depends on u and is wholly determined by that 
of V , Eq. (C77). It is interesting to note that, following 
Eq. (C108), M- 1 corresponds to Z\ | —1, and in the light 
of the considerations in Sec. C4a, the GS in the region 
|Z\| < 1 is disordered and the spectrum of the system is 
gapless. 

In the GS of the system corresponding to large 
values of g, ^{x, —ir) is locked into a minimum of 
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gcos(|3^(a;, — ir)) and the low-lying excitation spec- 
trum of the system is gapped. This spectrum is 
readily obtained by expanding to quadratic order the 
cos ((3^(a;, — ir)) in the action in Eq. (C90) around the 
above-mentioned minimum. The position xq of the min- 
imum of gcos(a;) depends on the sign of g; for .g < one 
has xq = (modulo 27r) and for g > 0, xq ~ tt (modulo 
2tt). Assuming that g > 0, on neglecting the —1 in the 
expansion of cos(a;) around x = n (by normal ordering) 
and effecting the transformation 



where 



S{x, -it) S{x, -it) + 
the action in Eq. (C90) reduces into 



S^S^+'—\g\l d^rf^(x,-zT), 



(C205) 



(C206) 



where we have used \g\ instead of g. This we have done 
on account of the fact that cos(a:;) = 1 — -I- 0{x'^) 
for a; — > (note the minus sign), whereby Eq. (C206) 
indeed applies also for g < 0. For completeness, the 
transformation in Eq. (C205) does not apply for the cases 
where g < 0. 

Expressing ^(.t, —it) in the double Fourier series 

f -*t) = ^ E e^^^^"-" kk, ^^n), (C207) 



where (cf. Eq. (C94)) 



iOn — — LUri 



(C208) 



with {ujn} Matsubara frequencies, and (cf. Eq. (33) and 
Eq. (C3)) 




(C212) 



is the mass gap aX k ~ ±kp. Since in this section we are 
tacitly assuming that > (see our pertinent remark 
following Eq. (C200) above), in what follows we shall 
denote \V\ hyV. 

In the case at hand, for the spectral gap £q at k = 
zb/cp one has 5g = 2M. As the physical spectral gap 
is small for Z\ i 1 (see Eqs. (22) and (C190)), for < 
Zi — 1 ^ 1 one can identify the u in Eq. (C212) with that 
specific to A = 1, that is it = nv-p/^ = nat (Eqs. (C126) 
and (C40)). Similarly for A', which for < Z\ - 1 <C 1 
takes a value satisfying < ^ — AT ^ 1, Eq. (C119); for 
a: < i the cosine term of the sine-Gordon Hamiltonian 
is relevant (see Eqs. (C95), (C97) and the subsequent 
remarks). Thus for < Z\ - K 1 from Eq. (C212) one 
obtains that 



M w 2V2t for < zi - K 1. 



(C213) 



Consequently, 



£,,{h) = 2M w 4V2i w 5.7 1. 



(C209) 



(C214) 

With 2t the bandwidth (for t' — 0), one observes that 
£q{L) is of the order of the bandwidth. 

We should emphasise that, in contrast to £G('*)e~'*, 
Eq. (C204), the quantity £q(1*) is physically meaning- 
less. This follows from the fact that in expanding the 
action as in Eq. (C206), we have taken no account of the 
reality that deep inside the massive phase (and indeed 
i is a very large mass), the parameter u, which 
is well-defined only in the region — 1 < Z\ < 1, carries no 
physically relevant information about the system under 
consideration. 

We are now in a position to consider the cases (1), (2) 
and (3) indicated above, on page 30. 



for the S in Eq. (C206) one obtains that 

S - ^^Yl[ojl+{ukf+\g\u^^'\S*{k,iu;n)S{k,t^r.), 

(C210) 



where 



i*{k,iujn) = i{~k,-iLOn)- 



The action in Eq. (C210) being quadratic in the modes 

{(p{k,iujn)}, the zeros of the expression enclosed by the 
square brackets in Eq. (C210), subsequent to the appli- 
cation of the analytic continuation 

w„ —iuj ± 0"'", y G R, 

coincide with the energy dispersions of these modes. One 
thus obtains that 



±{{ukf 



M 



2n1/2 



(C211) 



a. Case 1: zx ^ — zy 

In this case, one may identify the RHS of Eq. (C195) 
with zero (recall our convention, whereby z±_ > 0), im- 
plying that Z||(Z) is independent of I and thus equal to 
Z||. Using this result in the expression on the RHS of 
Eq. (C196), one immediately obtains that 



zO e^]' 



By defining li, as the value of I for which z±{l) 
more generally, z±{l) ~ 0(1)), we deduce that 

e-'*«|z"J-^/^% 
and thus, following Eqs. (C204) and (C212), 

£.^2A^«.4./^K|-/^". 
V Tra ' ' 



(C215) 
= 1 (or, 

(C216) 
(C217) 
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6. Case 2: z± = — zy 

111 this case, the differential equations in Eqs. (C195) 
and (C196) can be expressed as 



One readily verifies that 



Zq(0 



_L 



(C218) 



(C219) 



satisfies both the differential equation in Eq. (C218) and 
the appropriate boundary condition. Denoting, as above, 
the value of I for which z_l(Z) = 1 by i^,, from Eq. (C219) 
one obtains that 



1 = 



(C220) 



Hence, by Eqs. (C204) and (C212) one arrives at 

1 



2VKu 
= 27W = \\\ exp 



1 - 



(C221) 



where the last approximate expression is based on the 
consideration that in the case at hand z'[/|z|j| ^ 1 (re- 
call the above reasoning leading to Eq. (C222)). Conse- 
quently, following Eqs. (C204) and (C212), one obtains 
that 



£r, = 27W « 4 



2VKu 



exp 



y2(zO/|z||| 



1)^ 

(C226) 

This is the RG-based counterpart of the expression in 
Eq. (C190). The similarity between the two expres- 
sions clarifies our above choice for the condition (ij^/* = 
iTr^lzjjl. 

We note in passing that on identifying Zj^/|z||| with 
A = Vl(2t) for Z\ I 1, one can deduce the explicit ex- 
pression of the constant A, introduced in Eq. (C194), in 
terms of K\ for Z\ | 1, one obtains A - i7r^(l - IKf . 
Since for Z\ ], 1 one has K \ i, it follows that A vanishes 
in the limit of Zi = 1. 



c. Case 3: z± > — zy 

In this case (J.^ = V^M- > 0, which however close to 
transition it is nearly equal to zero. Consequently, on 
account of Eq. (C200), in the case at hand to a good ap- 
proximation one has z'J^/zjJ = — 1. Following Eq. (C203), 
one therefore to a good approximation has 

= TTi, (C222) 

as a result of which Eqs. (C201) and (C202) reduce into 



zii(0 = ^lCot(^l/), 



(C223) 
(C224) 



Inspired by the expression in Eq. (C190), here we define 
li, by the condition 

= ^Tr^lzjll = n\l - rf(0)/2) = n\l - 2K), 

where d(0) = = |3V(47r) = 4if, Eqs. (C197), (C95) 
and (C92). One thus has 



/ = _ = 'Jll ~ '1 

2y^(zO )2 - (z||)2 2y'2(zV|zlj|-l)' 

(C225) 



6. The single-particle Green functions Gj^-^{k;iuj„), 
A, Be {L,R} 

Here we determine the single-particle Green functions 
Gi^%{k; iuJn), A,B & {L, R}, as defined through Eqs. (31) 
and (32), in terms of the form-factors of the soliton- 
generating field operators pertaining to the quantum 
sine-Gordon Hamiltonian as determined by Lukyanov 
and Zamolodchikov (LZ) [32]. 



a. Preliminaries 



We first introduce (cf. Eqs. (C92) and (C77)) 



(p[x, -IT) 



^<P{x,-iT), (C227) 
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1 



aV 



(C228) 
(C229) 



whereby the action in Eq. (C90) can be written as 
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ydV{^[(9,^(x,-ir))' + (9,(p(x,-ir))'] -2/xcos(|3^(x,-»t))}. (C230) 

This Euchdcan action coincides with that underlying the considerations by LZ in Ref. [32]. The minus sign on the 
RHS of Eq. (C227) is not significant and serves only to bring our subsequent notation in conformity with that of LZ 
[32]. 

With reference to Eqs. (C82) and (C83), and in view of the fact that through the scaling transformations in 
Eq. (C88) and (C89) the K in the action in Eq. (C230) has been identified with unity, one has 

fl{x, -ir) ^ i dr^{x, -ir) = 9^(p(x, -ir). (C231) 

V87r 



Making use of this expression and that in Eq. (C227), taking into account the transformations in Eqs. (C88) and 
(C89), from the second expression in Eq. (C6) (in which we suppress the Klein factors C/± — see the remark in the 
paragraph following Eq. (027), on page 14), combined with Eq. (C27) (cf. Eq. (25)), one readily obtains that 



E{x, *r) 1 ^£^g^p[_ 1 r dx'9.cp(x',-*r)] exp[T^cp(.T + 0+,-*r)]. (C232) 



These expressions are to be compared with that of the 
non-local field 0^g{x) considered by LZ [32] and presented 
in Eq. (2.4) of Ref. [32] (here n is an index and is not to 
be confused with a power). From this comparison one 
deduces that in employing the form factors as calculated 
by LZ [32] in the context of our present considerations, 
one has the following relationships: 

L{x,-iT) J V27ra 4 

(C233) 

We note that in Ref. [32] x = (x, y), which is to be identi- 
fied with the Cartesian vector r = (x, t) in our notation. 
In the following we shall therefore denote the 0"{x) of 
LZ by d^ix,-iT). With reference to Eqs. (C85), (C86) 
and (C87), we point out that 

L{x, —ir) = L{x, —it), R{x, ^it) = R{x, —ir), 

(C234) 

so that, in view of Eqs. (31) and (32) and the second 
equality in Eq. (33) (for L 00), one has 

Gi,f^{k; iujn z) = -Gi ^(fc; -> - z), (C235) 
u - - u 

where w„ = w„/u (see Eqs. (C208), (C87) and (26)). 
The result in Eq. (C235) is a direct consequence of the 
Lorentz invariance of the continuum theory under consid- 
eration, so that it does not strictly apply for the CDW 
GS (see Eq. (C273) below and the subsequent remarks). 
Since the energy gap £c in the single-particle excitation 



spectrum of the system under consideration diminishes 
exponentially for V [ Vc, Eq. (22), for V sufficiently 
close to Vc the deviation of the function on the LHS of 
Eq. (C235) from its exact counterpart is however expo- 
nentially small. 

To determine the functional forms of the operators 
L{x,—iT) and R{x,~iT), Eq. (C234), one should reahse 
that for V ~ 0, R{x, —ir) is a function of z and L(x, —ir) 
a function of z, Eqs. (C46), (C50) and (C51); this remains 
the case for V ^ 0, and |Z\| < 1, barring the fact that 
for y ^ 0, and |Z\| < 1, the Wp in Eq. (C46) is trans- 
formed into the u in Eq. (C99). Similarly as regards 
^{Xj—ir) and 0{x,—iT), in the light of the expressions 
in Eq. (C47). Using the identities z = -^{t + ix) and 
z = ^iz ^ ix), one immediately infers that the expres- 
sion for i(x, —it) is deduced from that of L{x, —it) by 
substituting, in the explicit expression for the latter op- 
erator, x/u by X and r by r. Similarly for R{x,—iT). 
Conversely, the expression for L{x, —ir) is deduced from 
that of L(x, —it) by substituting, in the explicit expres- 
sion for the latter operator, x hy x/u and r by t. Simi- 
larly for R{x, -ir) (see Eq. (C270)). 

With \^nJ denoting the A^e-particle GS of Ti., we are 
to evaluate the following quantities: 

this on account of the fact that for r > one has (cf. 
Eq. (32)) 



-7ri(<;-?')/4 ^ ^ 

lim ^AB(a;,r) = _,v^^•^2t 

p^ca -'- 2'Ka 



-1, A,B = R, 



-1, A.B^L. 



(C236) 
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Following Eqs. (31) and (33), from these quantities one 
calculates G(fc; z) (see the following two paragraphs). Wc 
note that, without the '±' with which the expression on 
the RHS of Eq. (C232) (and similarly Eq. (25)) is pre- 
multiplied (see our remarks following Eq. (C27) above), 
the phase factor on the RHS of Eq. (C236) would have 
been e+"(^-'^')/4. 

The expression in Eq. (C236) is in general not suffi- 
cient to calculate the zero-temperature Green function 
G'a.b(2;, it^n), for the limit /3 — > oo in Eq. (31) is to be 
effected after evaluating the integral with respect to r. 
Similarly for Gx^Bi^] z). To be explicit, by effecting the 
limit /3 — > oo before evaluating the last-mentioned inte- 
gral, one can demonstrate that the resulting Green func- 
tion satisfies lm[Gi fi{k;Lu + ir])] = for all lo < when 
77 = 0^ . One can convince oneself of this fact by consider- 
ing Eq. (25.18) in Ref. [22], which reveals that evaluating 
the integral with respect to t over [0, 00) of the equiv- 
alent of the expression in Eq. (C236) is tantamount to 
identifying the (1 - n"j) in Eq. (25.18) of Ref. [22] with 
its zero-temperature limit, which is identically vanishing 
for all j corresponding to > /i. In contrast, use of 
the expression in Eq. (31) has the effect of totally com- 
pensating (1 — n°) by a factor —1/(1 — n^) which in the 
limit f3 — 00 is only operative for e'j > fi (in Ref. [22] , n'j 
denotes the non-interacting Fermi distribution function). 

In our case, the expression in Eq. (C236) contains suf- 
ficient information for calculating the zero-temperature 
Green function G^ Bi^^i^n), or Gi ^i^; z). This follows 
from the Lorentz invariance of the theory (which applies 
only approximately for V > Vc) whereby for A ^ B one 
has 

G'a,a(— fc; +^^) = +G||(fc;z), 

GA.li+k'j—z) = — G||(fc;z), 

Gi,|(-fc;-z) = -Gg,i(fc;z). (C237) 

These results are obtained from the spectral representa- 
tion of the Green functions Gi^^{k; z), A,B£ {L, R}, in 
terms of the single-particle spectral function A^j(fc;tj), 

Eq. (C287), and the symmetry properties (for A ^ B): 

Ai,A(k; (^) = A^^Bik; -w) = ABj(-fc; w), 

and 

^A,B(fc; ^) = ^B,A(-fc; -^^)- 

In what follows, we shall denote the functions 
Gf^Bi^T^^n) and GAj(fc;z) as deduced from the expres- 
sion in Eq. (C236) by G^. (x, iWn) and Gj^^{k]z) respec- 
tively. It can be shown that for A ^ B (see Eq. (C282) 
below) 

Gy-k;z) ^ G+g(fc;z), 

Gy-k;z) EE G+,(fc;z), 

G+g(-fc;z) ^ Gt^(fc;z). (C238) 
Making use the identities in Eq. (C238), one readily 
verifies that, for A 7^ B, the following Green functions 



satisfy the identities in Eq. (C237): 

G^j{k;z) = Gt^(fc;z)-G+g(-fc;-z), 
G:,.B{k;z) - Gf^{k; z) ~ Gl^{~k;-z). (C239) 

One can gain confidence that the Green functions as 
defined here are indeed the correct representations of 
Gp^x{k]z) and GiB{k]z) by identifying z with lj + irj^ 

Tj —> 0, and realising that Im[G^g(fc;[j -I- irj)] = for 

u} < and Im[Gt . (fc; -|- irj)] = Im[GA,B(^'; + *'?)] for 

uj>0, A,B e {L,R}. 

Wc note that, following Eq. (C237), from Eq. (33) one 
obtains that 

G(-/c;-z) = -G(fc;z). (C240) 

From Eqs. (C237) and (33) one further deduces that 

G(±fcp; -z) EE -G{±k^; z). (C241) 

This result is remarkable in that it establishes that for 
G{±kp; z) continuous at 2; = 0, one has 

G(±fcp; 0) = 0. (C242) 

This equality is to be compared with that in Eq. (9). One 
should also note the numerical results in Fig. 2. 

h. Technical details 

Let Z^{9) and Zl(0) denote the annihilation and cre- 
ation operators for the single-particle state character- 
ized by isotopic index e, e G and rapidity 9. 
These operators are not canonical, however satisfy the 
Zamolodchikov-Fadeev algebra [57] (also Ch. IV, Sec. 34 
in Ref. [4]). Here e = — refers to solitons (kinks) and 
e = + to anti-solitons (anti-kinks). Since \^n^) is a state 
free from solitons and anti-solitons, it is a vacuum state of 
the theory. Below we shall therefore employ the notation 

10) = lif-wj. (C243) 

One has 

Z,(6')|0) = 0. (C244) 

A multi-particle state, consisting of a total of n solitons 
and anti-solitons is determined according to 

10i,ei;...;0„,e„) = Zt^(^i)...Ztj0„)lO). (C245) 

The identity operator / of the Fock space under consid- 
eration can therefore be decomposed as follows [58]: 

X [01, ei; . . . ; 0„, e„)(0i, ei; . . . ; 0„, e„|. (C246) 
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We note that in the notation of LZ [32] the vacuum state |0) is denoted by | vac) . and a basis state containing n + N 
soHtons. with rapidities ^i,. . . , 0n+N, and N anti-sohtons, with rapidities O'l,. . . , by 

\A^i6^)...A^i9,,+N)A+i9[)...A+i6'^)). 

With E{9) and P{0) denoting respectively the energy and momentum of a sohton/antisohton corresponding to 
rapidity 8, one obtains that 



{WNMix,-lT)Ol\0,0)\^N, 



E E 

n = l {£!,...,£„} 



=° dOi... d0n 

oo (2^)"n! 
x(O|d2|0i,ei; 
I 



|2 



, ei ; . . . ; , e„ )*. (C247) 



For E{e) and P{e) one has (Eq. (2.5) in Ref. [57], and 
Eq. (34.1) in Ref. [4]) 

E{9) = m cosh(6i), P{9) = m smh{e), (C248) 

where m is the sohton/antisohton mass (see the re- 
mark following Eq. (35)). The expression in Eq. (C247) 
makes explicit the way in which the Green functions 
^Ab^^^' ^'i^")' 4, ^ S {L,U}, are determined in terms 
oi form factors (OlO^l^i, ei; . . . ; 0„, e„), s& {-|,+|-}. 
Conservation of topological charge implies that in the 
case at hand only those form factors are non-vanishing 
for which number of solitons is by two units in excess of 
number of anti-solitons; in general, only those form fac- 
tors corresponding to Of are in principle non-vanishing 



for which number of solitons is by p units in excess of 
number of anti-solitons. Consequently, the first non- 
vanishing term on the RHS of Eq. (C247) corresponds 
to n = 2. By the same reasoning, we have already sup- 
pressed the contribution of |0)(0| in / to the expectation 
value in Eq. (C247). 

Close to insulator-metal transition (i.e. for A | 1), 
the first sum on the RHS of Eq. (C247) is dominated 
by the terms corresponding to n = 2. By the reason- 
ing of the last paragraph, in this region only the exited 
states corresponding to two solitons, i.e. to ei = — and 
62 — contribute to the sum Y^r i on the RHS of 
Eq. (C247). For Z\ J, 1 we therefore employ the following 
approximation: 



One has (Eq. (2.12) in Ref. [32]) 



1 clf) clO r 1 

J ^^^cxp[^^(^P(0,)^-m)T)](O|af|0l,-;^^2,-)(O|a,?,|0l,-;02,-)* 



(C249) 



|a2|0i,-;(?2,-) = v/^exp[^] exp[-i(0i + 02)]c;(0i-02) 



(3 



= ,/z^)cxp['-^] exp[^{e, +62)] g (9,-92), s = ^^, {-1,+1}, (C250) 



where Z2{s) = Z2{<,) is a normalization constant, and (Eq. (2.13) in Ref. [32] — see appendix D) 

At sinh^((l - i9/T:)t) sinh((£, - l)t) 



g{9) ^ iCi sinh(6l/2) cxp 
in which (Eq. (2.14) of Ref. [32]) 

Ci = exp 



sinh(2t) cosh(t) sinh(y) 



dt sinh^(i/2)sinh((l - l)t) 



t sinh(2t) cosh(t) sinh(£,t) 
I 



(C251) 



(C252) 



and (Eq. (2.15) in Ref. [32]) where we have used Eq. (C228) above. The quantity £, 

is not to be confused with the correlation length ^. The 

' " -^2 ^ (C253) 



1 - 1-2K 
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coefficient Ci in Eq. (C251) is a function of £,. Since 
K t i for Z\ I 1 (Eq. (C97) and the surrounding text, as 
well as Sec. C4c), it follows that 



£, T +00 for Ai 1. 



(C254) 



Since our interest is restricted to the region where J, 1, 
in what follows we shall employ the following simplifying 



result: 



sinh((t-l)t) 
lim . 

£,-►00 smh(£,t) 



t > 0. 



(C255) 



For t > and £, finite, the leading-order correction to the 
RHS of this expression is equal to — e~^^*. 



J 



LZ [32] have conjectured an expression for Z„(s) which has been presented in Eq. (3.1) of Rcf. [32]. For the specific 
case of 71 = 2, one has (appendix E) 



0Fmr(| + f) 



d{s,2) 



exp 



\lt ( 1 / cosh(£,t) c-^(^+^)t-i X 
t l4sinh(^t) Vsinh((^+l)i)cosh(t) J~ ^ ° 



(C256) 



where (Eq. (2.8) in Rcf. [32]) 

!2 



2P^ 



One has 



8 2(3 



(3 5 

(i(± — , 2) ~ - as £, — > 00. 



^2 for s^±j- (C257) 



(C258) 



From Eqs. (C256) and (C257) one observes that for the 
case at hand, Z^i^s) takes the same value for s = +P/4 



and s 



-p/4. In appendix E wc consider the behaviour 



of Z2{s) and in particular show that for a constant 
value of the soliton mass m, this function decays like 
as 00. 

As is evident from the expressions in Eqs. (C249) and 
(C250), for our considerations we need to calculate |tJ(0)|, 
and not G{0). Consequently, we can simplify the calcula- 
tions to be performed by separating the pure phase factor 
that may contribute to Q{9). To this end, we employ the 
identity 



sinh^((l - i6/TT)t) = - siw^ [Ot / -k) cosh(2i) +sinh^(t) - i sin(6't/7r) cos(6't/7r) sinh(2t). 



(C259) 



One observes that for 6', t £ M, the last term on the RHS of this expression gives rise to a pure phase in the expression 
for Q{9). Since the contribution to Q(9) of the second term on the RHS of Eq. (C259) is independent of 9, we absorb 
the contribution to Q{9) as arising from this term into Ci and thus write 



where 



Ci = Ci exp 



g{9) = iCiZ^(6l)sinh(6l/2)cxp 
~ iCiU{9)smh{9/2)cyij> 

""dt sinh^(t)sinh((£,- 



At sin2(6't/7r)sinh((£,- l)i) 



t tanh(2i) cosli(t) sinh(£,<) 
°° dt sin2(6't/7r) c"* 





t tanh(2f) cosh(i) 



for £, 00, 



Q t sinh(2i) cosh(t) sinh(y) 



06 



21/6 



Ci = 1.356133165 ••• X Ci for 



(C260) 



(C261) 



U{9) = exp 



'dt sin(6't/7r)cos(6l</7r)sinh((£,- l)t) 



t 



cosh(i) sinh(£,t) 



(C262) 



in which © = cxp(l/12 - C'(-l)) = 1.282 427129 ... is the Glaishcr number; here C'(^) is the derivative of the 
Riemann zeta function (Ch. 13 in Ref. [37]). Clearly, \U{9)\ = 1, V0 G R. For £, ^ 00 one has 



Ci ^ exp 



'dt sinh^(V2) e" 



g t sinh(2t) cosli(f;) 



= 1/4 



21/6(33 Y pi^l 



^' = 0.932 937401 



(C263) 
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Following Eq. (C261) and (C263), one thus has Ci - 1.265 187 350 ... for f, c». 
Making use of the expressions in Eqs. (C249) and (C250), one obtains that 



{1'N^o'i{x,~lT)^l) {oM'^n:} ~ As,s' I d(?id(?2 \g{ei - 02)1 



X exp 



m(i[sinh(6'i) + sinh(6'2)]a; - [cosh(6'i) + cosh(6'2)]T) + 



where 



-4, 



\/Z2{s)Z2{s') rTTi . ^2(5) „ 



+ 1, s = s', 

±i, s = — s' = ±-x. 



(C264) 



(C265) 



In writing the last expression, we have made explicit the fact that ^JZ^Xs) takes the same values for s = +(3/4 and 
s = —^jA: (see our remark subsequent to Eq. (C258) above). In appendix D we present the closed expression for 
\Q{&)\ specific to large values of £,. 

Introducing the auxiliary variable B = Oi — 62, and using the transformation cxp(0i) — ^ y, the expression in 
Eq. (C264) can be expressed as (with s = <; |3/4, s' = ^' P/4) 



{'FNjo-iix,^iT)dl\o,o)\'FNj « As,s' I de \g{9)\-' i'^>{e-x-iT) 



1 2 -7-(m) , 



where 



4"j(0;a;,-*r) = e'^^+^'^^Z^ /"dy exp[mfz[^^^ + 

7n I \ ^ 2y 



, .r ^^-l , y^exp(-2^)-l ^ .2/2^1 , y^eM-'^0) + \^ 

tn|l| 1 ; J2^— |_- 



2yexp(- 



2j/ 



2yexp(- 



-(^ 



i^(?+^')/2(2mcosh(6l/2)Vr2 + a;2), r > 0, 
I 



(C266) 



(C267) 



in which K^{z) is the modified Bessel function of the lows that the difference between the X^"^) (6*; x, —it) cor- 
responding to = — 1 and that corresponding to = +1 
is entirely due to the function pre-multiplying the Bessel 
function on the RHS of Eq. (C267). From the last ex- 



second kind (Ch. 9 in Rcf. [23]). With 



(C268) 



pression in Eq. (C267) one trivially deduces that 



on account of ^, E { — 1, +1}, one has 
{-1,0,+1}. 



{9; X, ~^T) ^ {6; x/u, ~zr) = X^^' {9; x, -it). 

(C269) (C270) 
By introducing the function (see Eqs. (31), (33), (C94) 
Since K-,,{z) = K^{z) (item 9.6.2 in Ref. [23]), it fol- above as weh as Eq. (25.14) in Ref. [22]) 



I ,{9;k,iujn) = lim / dr / 

00 /n J-L/2 



from Eqs. (C236) and (C266) one obtains that 



dx e*(y"T-M 4!^)(0;a;,-^r), 



(C271) 



d^ |g(^^)|2iJ!(0;A;,iy„), ^,<r' 



-1, A,B^R, 
+ 1, A,B = L, 



(C272) 



where Z2 is the short-hand notation for ^2(5), which, as we have indicated following Eq. (C265) above, takes the 
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same value for s = +|3/4 and s = — 13/4. In arriving at 
Eq. (C272), we have made use of the fact that \G{0)\ and 

^/ (0; fc, iujn) are both even functions of 9. 
' Since w„r = w„t, Eqs. (C87), (C208), from Eq. (C270) 
one deduces that 

i^7^(0;fc,ic^„) EE -l['l]{9;k,iuj,,), (C273) 

from which and the expressions in Eq. (C239) one obtains 
the equality in Eq. (C235). The resuh in Eq. (C273) is 
the manifestation of the Lorentz invariancc of the con- 
tinuum theory under consideration. Although this in- 
variancc breaks down in the CDW GS, where m 7^ 0, we 
have not taken explicit account of this breakdown by the 
consideration that, for 1/ > and V sufficiently close 
14, m is exponentially small, a fact that we have empha- 
sized in our remarks preceding Eqs. (30) and (31). As 
regards the relationships between 



■,k,iuJn) and i[J,{9;k,iujn) 



with respectively 

Gi^^(k;iujn) and G-^^^{k;iuJn), 

we refer the reader to the remarks in the paragraph fol- 
lowing Eq. (C235), on page 33, which underlie the defin- 
ing expression in Eq. (C270). 

We evaluate the double-integral in Eq. (C271) by em- 
ploying the circular coordinates {g,ip)' 

TT 7T 

T^gcos{ip), x^gsm{ip), [--,-], (C274) 

whereby from Eq. (C271) one obtains that 



where 



-4 d r'^ 



t/2 



d^ c''^^/(^)(0;fc,iw„;^) 



(C275) 



l + (e)' + e(f -sinh-\f)) 



A + (f)2 



(C276) 



f Jl + (|)2 + f +zsinh-i(|) 
r±{{9]k,iujn;v) = , , 

(C277) 

a = a{iujn, k; ip) = ujn cos((p) — k sin((^), (C278) 



X = x{0) = 2m cosh(e'/2). 



(C279) 



In Eqs. (C276) and (C277), sinh~^(z) stands for the 
inverse of sinh(z), for which one has (item 4.6.20 in 
Rcf. [23]) 



sinh ^(z) = ln(z + Vl + z^). 



(C280) 



The property f[^l 



f'^i is a consequence of K-^{z) 



K^,^{z), to which we referred following Eq. (C269) above. 

From Eqs. (C275), (C276) and (C277) one trivially de- 
duces that for 1/ E { — 1,0,+!}, and any finite value of 
\iuJn\, one has 



J^^!(6';fc,za;„) 



as 



1^1 



(C281) 



In view of the leading-order asymptotic behaviour of 
\G{9)\ corresponding to \9\ 00 (see Eq. (D19) below), 
one infers that for |0| — > cxd the integrand of the integral 
on the RHS of Eq. (C266) decays like \9\'^/^ e^l^l/^. 

The order of differentiation with respect to A and in- 
tegration with respect to ip in Eq. (C275) is significant; 
exchanging this order, the integral with respect to (p may 
not exist. In our calculations, we have evaluated both 
the integral with respect to ip and the derivative with 
respect to A in Eq. (C275) numerically. For the latter, 
we have employed the central finite difference method as 
presented in item 25.3.21 of Ref. [23] (using 2h < 10"'', 
for which reason we have evaluated the integral with re- 
spect to V? to a relative accuracy better than (2/i)^). In 
Fig. 2 we present some numerical results of our calcula- 
tions based on the details presented in this section. 

In closing this section, we remark that through the 
variable transformation p —p in Eq. (C275), one im- 
mediately deduces that, with f == — 

1^^^ {9;k,iujn) = l''^^{9;~k,iujn), 

i'^^^{9;k,iujn) = i'^^^{9;-k,iLJn), 

j["\0;fc,zy„) EE xJ;^-^(0;-fc,zy„), (C282) 

which; following Eq. (C272), lead to the symmetry rela- 
tions presented in Eq. (C238). 



c. Some general remarks 

For |z| — *■ 00 and | argz| < ^tt, one has the asymptotic 
expression (item 9.7.2 in Ref. [23]) 

^ 7-1 , (7-l)(7-9) , 1 



2z I 



1 



2z " L ■ 8z 2!(8z)2 

(C283) 

where 7 ^ Av^. This asymptotic expression reveals that 
the infinite series in powers of l/(8z), enclosed by curly 
brackets, is terminating for 



7 = 1 



2 



(in which case it only consists of the single term 1), 

7 = 9 < 



3 



(in which case it consists of two terms, 1 and (7— l)/(8z)), 
etc. In our case, where v takes the tree values indicated 
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in Eq. (C269), the function i^(^_|_^/)/2(^) = K^{z) (as 
encountered on the RHS of Eq. (C267)) is not expressible 
as a product of z~^/^ e~^ and a finite series in powers of 
l/z. 

If one were to evaluate Gi^^{k\ z) for small values of \k\ 
and jzl, it would be reasonable to approximate K^{z) by 
the leading-order contribution on the RHS of Eq. (C283), 
an approximation whose validity is conditioned on 

2mcosh(6i/2)v/r2 + 2)2 > a. (C284) 

Employing this approximation, and using the identity 
(which underlies the expression 2.580.1 in Ref. [59]): 

a + bcos{ip) + c sin((y9) = (a + p) - 2p sin^ [ip) , (C285) 

where 

ip = 2i' + i}, taii(d) = ^, p = bsec(i}), (C286) 



one can express ^!(0; fc, iuJn) fully analytically in terms 
of the elliptic integrals of the first and second kind, 
F{iP,k) and E{ip,K) respectively (items 8.111.2 and 
8.111.3 in Ref. [59]). We note in passing that although 
the real part of = tan~-'^(c/6) is determined up to the 
constant mr, where n = 0, ±1, . . . , this ambiguity in the 
value of is of no consequence to the validity of the iden- 
tity in Eq. (C285). We do not present this closed expres- 

sion for ^i(6]k,iujn) here essentially for the fact that 
evaluation of elliptic integrals for complex arguments and 
moduli is not computationally less time-consuming than 
evaluation of the integral with respect to ip in Eq. (C275) 

(additionally, the full closed expression forX^ ^i{0\ k, iuJn) 
is both very extensive and involves some detailed pre- 
scriptions for singling out the appropriate branches of 
the many-valued functions -F(V', i^) and E{il)^ k)). 

It should be noted that since tn diminishes exponen- 
tially for V I Vc{t,t'), Eqs. (22) and (35), the validity 
of the condition in Eq. (C284) for all 9 implies that as 
V i Vc{t,t'), the approximation for Gi^^{k;z) based on 
the leading-order term in Eq. (C283) is accurate in princi- 
ple only for exponentially small neighbourhoods of A: = 
and z = 0. 



d. The leading- order term in the asymptotic series 
expansion of G{k; z) for \z\ — > oo 

The exact leading-order term in the asymptotic series 
expansion of the exact G{k] z) is l/z (in the units where 
ft = 1), where the coefficient 1 arises as a consequence 
of the following three properties (Sec. B.3 in Ref. [7]): 
(i) the completeness of the intermediate states used in 
the calculation of G(fc; z) (compare with the decomposi- 
tion of the identity operator / in Eq. (C246)), (ii) the 
anti-commutation relation [ak,a\]- = 1 for canonical 
fermion operators, Eq. (C4), and (iii) the normalisation 



{^N^\^N^) = 1- Equivalently, with 

A(fc;tj) = =F-Im[G(fc;w±iO+)], a; e R, (C287) 

the single-particle spectral function, the 1 in the above- 
indicated l/z is exactly the same constant as the 1 on 
the RHS of the following exact sum rule: 




dw A(fc;w) = 1, Vfc e IBZ. (C288) 



It follows that within the framework of our formalism, 
where we have replaced the sum with respect of n on the 
RHS of Eq. (C247) by the underlying summand corre- 
sponding to ri = 2, the coefficient of the leading-order 
asymptotic term of G(fc; z) corresponding to ]z] ^ oo is 
less than unity. To appreciate this fact, since A{k] lo) >Q 
for all Ld (a consequence of the assumed stability of the 
GS of the system under consideration), restriction of the 
last-mentioned infinite sum to its summand correspond- 
ing to n = 2, amounts to the imposition of a specific 
(soft) cut-off on A{k]Ld), which in general leads to the 
integral in Eq. (C288) taking a value less than unity. 

Due to the complexity of the full analysis of the asymp- 
totic behaviour of the Green function G(k\ z), for \z\ — > 
oo, as calculated here, Eq. (33), below we only present 
an outline of the full analysis, which we do not claim 
to be rigorous. The aim of this investigation is to show 
that the Green function G(fc; z) as calculated in this pa- 
per indeed decays to leading order like G/z for \z\ oo, 
where, following our above considerations, G is a real con- 
stant satisfying < G < 1 (the crucial property G > 
should be noted). To this end, below we consider the 
leading-order term in the asymptotic series expansion of 

2^ ^,{0; k,iujn) corresponding to \iuJn\ —>■ oo. In this con- 
nection, note that since l/z is analytic everywhere in the 
complex z plane outside z = 0, the direction along which 
we let z approach the point of infinity of the z plane is of 
no significance to our present leading-order calculations. 

For ip 7^ T7''/2 one readily verifies that dfj)'^^ /dX de- 
cays to leading-order like l/{iujn)'^ as \iuJn\ — ^ oo. It 
follows that the leading-order term of the asymptotic se- 
ries expansion of for \iujn\ ^ oo originates from 
the contributions to the integral with respect to tp in 
Eq. (C275) as arising from the small neighbourhoods of 
q=7r/2 inside [— 7r/2, tt/2]. This is directly plausible, since, 
following the first expression in Eq. (C274), (p T7i'/2 
corresponds to r — > for any finite value of g. Expanding 
cos{p) and sin((/3) around p = ^t:/2, from the expression 
in Eq. (C278) one deduces that 

--^(±^+;)±- for ^^t|, (C289) 

from which one infers the following leading-order asymp- 
totic expression for \iujn \ — *■ oo: 



40 



-4 d r 



7r/2+x/|iy„ 
-7r/2 



+ C2 



/2 

7r/2-x/|jy„| 



(C290) 



r 



where Ci and C2 are constants of the order of unity. 
Since | fl^^ \ is bounded and non- vanishing for ip over the 
intervals of integration indicated, and since the lengths 
of these intervals is equal to x/|zu;„|, one immediately 
observes that indeed the RHS of Eq. (C290) decays like 

~ (m) 

x/{i^n) for \i^n\ ~^ 00. This in turn implies that ^, 
to leading order decays like l/(«WnX) for \iujn\ 00. We 
have verified the correctness of this result numerically. 

With reference to Eq. (C281), one observes that the 
processes of effecting the two limits |0| — > 00 and \iujn\ 
00 do not commute: our considerations in the previous 
paragraph have revealed that for any finite value of 9 (and 
therefore of x), the in Eq. (C281) is replaced by 

l/(iw„x) K^nl approaches 00. This ambiguity is by- 
passed through the observation that for any finite I^Wnl, 
the integrand of the integral in Eq. (C272) decays like 
|0|i/2g-|ei/2 fQj. ^ gqg^ (C281) and (D19), so 

that to exponential accuracy the upper bound of the in- 
tegral on the RHS of Eq. (C272) can be replaced by a 
finite (large) constant 9q (in this connection, one should 
note that, following Eq. (C279), x increases exponen- 
tially with as 1^1 — > 00). With reference the conclusion 
arrived at in the previous paragraph, on replacing the up- 
per bound of the integral on the RHS of Eq. (C272) by 
a finite constant, the conclusion that for \iujn\ — > 00 the 
corresponding ^{k]iujn) decays like 1/(1^,1) is imme- 

I 



diate. With reference to Eq. (C239), we thus conclude 
that for \z\ 00 the corresponding G{k;z), Eq. (33), 
indeed to leading order decays like C/z, where C > 0. 
Our numerically-calculated G{k; uj±iO'^) (see Fig. 2) con- 
form with this observation; explicitly, our numerical cal- 
culations show that for uji < \u;\ < UJ2 this function is 
accurately described by a function decaying like I/lo"', 
where a increases monotonically towards 1 for lui <C uj2 
as wi, ti;2 — *■ cxD. 



APPENDIX D: THE CLOSED FORM OF \g{e)\ 

FOR f, ^ 00 

The leading-order expression for G{0) in the asymp- 
totic region corresponding to f, — > 00 is presented in 
Eq. (C260). Here we determine the closed form of the 
absolute value of this expression. To this end, we define 



/ dt sin2(0V^) e-* 



t tanh(2t) cosh(t) ■ 



(Dl) 



Making use of the expressions for tanh(2t) and cosh(t) 
in terms of y = e^* and y~^, followed by a partial-fraction 
expansion of the resulting expression, one deduces the 
following identities: 



tanh(2i) cosh(t) 



1 1 



2y-l 



1 1 



2 y+l 



2 VI 



1 



1 



4sinh(i) 4cosh(t) 1 + cosh(2i) 2sinh(2i)' 
I 



(D2) 



One has 



Further, 



dt siT?{et/TT) c- 
t sinh(i) 



In 



sinh(0) \ 



(D3) 



where we have used the analytic continuation of expres- 
sion 3.553.1 in Ref. [59] for a < 1 to a along the imagi- 
nary axis (explicitly, to a = i^/tt, where G R). From 
Eq. (D3) one immediately obtains that 



dt sv[?{et/T:) c-2* 
~t sinh(2t) 



In 



/sinh(6i/2)^ 



61/2 



(D4) 



dt sin^(e't/7r) e"' 
t cosh(i) 



In 



- coth ( - 
2 4^ 



(D5) 



where we have used the analytic continuation of expres- 
sion 3.555.2 in Ref. [59] for a < 5 to a along the imaginary 
axis. 

In order to obtain the closed expression for the integral 
corresponding to the third term on the RHS of Eq. (D2), 
we introduce the following auxiliary integral (this would 
not have been necessary, had the expression 3.557.1 in 
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Rcf. [59] been correct): 



/o t 1 + cosh(t) 
One has 

dt sm^{et/TT) 



(D6) 



---An3+T{--^)]. (D7) 



t 1 + cosh(2i) 4 V ^t: 



For determining the closed expression for T(p), it 
proves convenient to consider dT(jj) /dp and subsequently 
obtain T{p) through integrating the latter function. One 
obtains (see expression 3.541.8 in Ref. [59]) 



dp 



—Tip) = dt —J- = 1 - 2pm, (D8) 

cosh (t) 



where (Sec. 8.37 in Ref. [59]) 

1 / , + ^ 



P(^)=2(^(- 



(D9) 



in which ipi^) is the digamma function (Sec. 8.36 in 
Ref. [59] and Sec. 6.3 in Ref. [23]) 



d 



V'(z)--ln(r(z)). 



(DIO) 



In the following we shall encounter the polygamma func- 
tion ■!/''^™)(z) defined according to (Sec. 10.2 in Ref. [60]) 



(Dll) 



Thus V'(z) = i^^°'>{z). One has 



V-('"-i)(z)= / dz'V^"^(^')- (D12) 
Jo 

For completeness, we point out that the digamma and 
polygamma functions in Ref. [60] are defined in terms of 

r{z + i) = zr{z), 

whereby the polygamma function Fm{z) in Ref. [60] is 
identical to ip^'^^z + 1). Further, although ip'^^^^z) is 
defined for integer values of m, this function can be ana- 
lytically continued into the complex ly plane whereby one 
can deal with ■0'^'')(z), where i/ e C [61]. 

Since, following Eq. (D6), T(0) = 0, on integrating 
both sides of Eq. (D8) with respect to p over [0,p], one 
obtains that 



Tip) = p - 2p(ln (r(^)) - In (r(|))) + (^) - ^/-'^ (|)) - In {2^/\&^), 



p+l. 



(D13) 



where © = e^/^^^^'^-^) = 1.282 427 129 .. . is the Glaisher number. In arriving at the expression in Eq. (D13), we have 
determined the integral dp' p'P){p') through applying integration by parts, making use of the defining expressions 
in Eqs. (D9) and (DIO), and of the exact result 

''dy {ln(r(^)) -ln(r(^))} =2(^(-)(^) -^(-)(|)) - iln(25/W). (D14) 

Collecting the closed forms of the relevant integrals of the terms on the RHS of Eq. (D2), one arrives at 



^ l^^/ 8sinh(0)sinh^(g/2) x 1^^,10 lO 



614 cosh(6'/2) / 4 



One thus has 

g[e) ^ iCiU{e)suMd/'^) e^-^^^^ ^iCiU{9)(^- 

Since 



'sinh^(6i/2)cosh(6'/2)\i/8 



L[r{-)+T{-- 



8 sinh(6') 



ln(2) 
27r2 



exp 



(n- 



T{ 



(D15) 



for £, — !■ oo. 
(D16) 



0, 



l+ln(^) 



(D17) 



■ln(6'2) + — — 6'^±oo, 



42 



with 

1 



sinh^(6'/2)cosh(6'/2)\i/8 J 2 ' ' 



8 sinh(6l) ) 



(D18) 



it follows that 



l\Ci\\e\ 



{en) 



1/4 , 



-\c,\\e\ 



1/4 p|e|/4 



0, 



6* ^ ±oo, 



(D19) 



22/303 

where, following Eq. (C261) and (C263), one has 



Ci - 1.265187 350 



for £, —t oo. 



(D20) 



Recall that \U{e)\ = 1 for 6* G M. The cusp in \g{e)\ 
at = is entirely due to the function in Eq. (D18), 
so that it is also present in Q{9) at 6 = 0. Explicit 
calculation shows that for 16*1 < 2 {\6\ > 2)the upper 
(lower) leading-order expression in Eq. (D19) relatively 
accurately describes the behaviour of the exact |^(6')|. 

For completeness, we point out that in our numerical 
calculations we have employed the following approxima- 
tion for i/j^"^) (z): 

« ^,,(1 - |z|) /i(z) + i3„{\z\ - 1) /2(z), (D21) 



where 



exp(— x/?;) + 1 ' 



(D22) 



and fi{z) and /2(z) arc the asymptotic series expansions 
of V'^"^'' (-2^) corresponding to z — > and z ^ oo respec- 
tively. The Fermi function in which < 7y ^ 1, 
serves to prevent discontinuity in the function on the 
RHS of Eq. (D21) on the circle \z\ = 1 (in practice, 
we have used rj = 10^"^). For our calculations, we have 
equated fi{z) and /2(z) with the relevant fifteenth-order 
asymptotic series expansions of ^^~^\z), which we do 
not reproduce here; these can be readily obtained with 
the aid of Mathematica. We note that both scries ex- 
pansions involve ln(z) so that they are not Taylor expan- 
sions around z = and z = oo. We further note that 
the approximation in Eq. (D21) is suited for our present 
calculations (where 6* S R) and should be used, if at all, 
with care in other applications. 



APPENDIX E: ON THE BEHAVIOUR OF a/I^sJ 

As we indicated in appendix C, in Ref. [32] LZ have 
conjectured an expression for ^JZ^J^ (see Eq. (3.1) 
in Ref. [32]). We have reproduced this expression in 



3 - 



N 
>• 




200 



FIG. 3: (Colour online) The log-log plot of ^/Z^ for 
s = ±|3/4 and m = 1. Solid line represents the numeri- 
cally calculated \J Z^is) and the broken line the expression 
yl/£,^/^^, where A = 4.947 289 ... (we have deduced both 
this value and the power 3/16 from our numerically-calculated 
^ Z2{s); numerically, we reproduce the exponent 3/16 to an 
accuracy of 1 part in 10^). For arbitrary m, the corresponding 
\fZj2(s) is obtained by multiplying the data presented here 



by m' 



d(s,2) 



where d(s, 2) 



Eq. (C257), in which P = 2^" 
has d(±P/4,2) ~ 5/8, Eq. (C258). 



+ 1/(2^=) for s 
£./(£. + 1). For f. 



±P/4, 



Eq. (C256) for the specific case of n = 2 and s = ±|3/4. 
In Fig. 3 we display this function for the specific case of 
m = 1 as calculated numerically; in the underlying cal- 
culation, we have taken full account of the dependence 
of C\. Eq. (C252), on £,. The value for m that is relevant 
for the case of Z\ | 1 is deduced from the expressions 
in Eqs. (22) and (35) (sec Eqs. (C190) and (C226)). In 
Eq. (35), u is the renormalizcd Fermi velocity at fc fcp 
in the limit A\\\ according to Eq. (0123), in this limit 
one has u=\v^ (Eqs. (C40) and (C126)). In this con- 
nection, we point out that in Ref. [62] (Eq. (12) herein) 



2r(f) 




l)/2 



(El) 



where is given in Eq. (C229); we have replaced the 
^l in the original expression by since sign of /x can 
be altered by means of the transformation in Eq. (C80). 
From the expression in Eq. (El) one readily obtains that 
(items 6.1.33 and 6.1.40 in Ref. [23]) 



m ■ 



2V2e 

eT 



£./2 



as f, — > oo. 



(E2) 



43 



where 7 = 0.577 215 ... is the Euler constant. This 
stronger than factorial growth in m for f, — > 00 is what 
one expects in the strong-coupHng (SC) regime ^ — Z|| 
(Sec. C5a). To appreciate this fact, one has to consider 
the RG result in Eq. (C217) (recall that m = M/{2u)) 
and take into account that from Eqs. (C198), (C95), 
(C228) and (C253) one has 

_ 1 = 1 = 1 ^ 1 ^ £.+ 1 

zjj " 2- d" 2(1-132/(8^)) " 2(1- " 2 

(E3) 

With reference to Eq. (C217), note that (Eq. (C194)) 

z'[ = Zi(/ = 0) = = 0), 

where z{l = 0) = z is the dimensionless coupling constant 
presented in Eq. (C191); one observes that indeed z^ is 



indeed proportional to |/.t|, Eq. (C229). 

From Fig. 3 it is evident that for a constant m (i.e. 
independent of £,) y^Zj^ipT^) decreases like l/i^/"^^ for 
£, 00; more precisely, for £, > 3, \/ Z2(±P/4)|m=i turns 
out to be relatively accurately described by the power low 
A/f,3/i6^ with A = 4.947289 .... Our numerical calcula- 
tions reveal that for a constant m, \/ Zi(±|3/2) increases 
like £,3/16 for — > 00. Interestingly, from the explicit 
expression in Eq. (4.24) of Ref. [32], one immediately ob- 
tains that, for m constant, ^4(0) decays like l/£, for 
£, ^ (X). It follows that unless m is made to vary with 
an appropriate power of the conjectured expression for 
the normalization constant ^J Zn{s) is cither unbounded 
or vanishing in the limit £, = 00; in both cases, £, = 00 
amounts to a singular point of the theory. □ 
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